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Abstract 

This paper explores and ties together three themes. The first is to establish regular¬ 
ity of a metric tensor, on a manifold with boundary, on which there are given Ricci 
curvature bounds, on the manifold and its boundary, and a Lipschitz bound on the 
mean curvature of the boundary. The second is to establish geometric convergence 
of a (sub)sequence of manifolds with boundary with such geometrical bounds and 
also an upper bound on the diameter and a lower bound on injectivity and bound¬ 
ary injectivity radius, making use of the hrst part. The third theme involves the 
uniqueness and conditional stability of an inverse problem proposed by Gel’fand, 
making essential use of the results of the hrst two parts. 


1. Introduction 

The goals of this paper are to establish regularity, up to the boundary, of the met¬ 
ric tensor of a Riemannian manifold with boundary, under Ricci curvature bounds 
and control of the boundary’s mean curvature; to apply this to results on Gromov 
compactness and geometric convergence in the category of manifolds with bound¬ 
ary; and then to apply these results to the study of an inverse boundary spectral 
problem introduced by I. Gel’fand. 

Regularity of the metric tensor away from the boundary has been studied and 
used in a number of papers, starting with [DTK]. The tack has been to construct 
local harmonic coordinates and use the fact that, in such harmonic coordinates, the 
Ricci tensor has the form 

(1.0.1) ^s') 2 Ric^ni ■ 

Here A is the Laplace-Beltrami operator, applied componentwise to the components 
of the metric tensor, and is a quadratic form in with coefficients that 

are smooth functions of Qij as long as the metric tensor satishes a bound Ci\ri\^ < 
gjk{x)r]^r]'^ < C2\r]\‘^, with 0 < Ci < C2 < 00. If one is given information on the 
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Ricci tensor, one can regard (1.0.1) as an elliptic PDE for the metric tensor, and 
obtain information on its components, in harmonic coordinates. 

The notion of compactness of a family of Riemannian manifolds and of geometric 
convergence issnes from work of J. Cheeger [Ch] and M. Gromov (cf. [Gr], the revised 
and translated version of his 1981 work). The role of harmonic coordinates in the 
stndy of snch geometric convergence has been exploited in a nnmber of papers. It 
was nsed in [P] and in [GW] to obtain a compactness resnlt, assnming a bonnd on 
the Riemann tensor, and some other geometric qnantities. In [Anl] there was a 
snccessfnl treatment of compactness given a snp norm bonnd on the Ricci tensor, 
an npper bonnd on the diameter, and a lower bonnd on the injectivity radins, for 
a family of compact Riemannian manifolds of a hxed dimension. Gonvergence was 
shown to hold, for a snbseqnence, in the G^-topology, for any r < 2. (A dehnition 
of geometric convergence is recalled in §3.) 

One of onr motivations to extend the scope of these resnlts to the category of 
manifolds with bonndary arises natnrally in the stndy of a class of inverse problems. 
In these problems, one wants to determine the coefficients of some partial differ¬ 
ential eqnation in a bonnded region via measnrements of solntions to the PDE at 
the bonndary. Snch problems arise in varions areas, inclnding geophysics, medical 
imaging, and nondestrnctive testing. One problem, formnlated by I. Gel’fand [Ge], 
consists of Ending the shape of a compact manifold M with bonndary dM and the 
metric tensor on it from the spectral data on dM. Namely, if R\ is the resolvent of 
the Nenmann Laplacian on M, the Gel’fand data consists of the restriction of 
the integral kernel R\{x,y) of the resolvent io x,y G dM, as A varies over the re¬ 
solvent set of A-^. Another formnlation of Gel’fand’s inverse problem will be given 
in §4. 

For snch an inverse problem, the first issne to investigate is nniqneness. In the 
context of metric tensors, this was established for the Gel’fand problem in 
[BKl] taking into acconnt the nniqne continnation in [Taj. See also [Bz], [NSU], 
[Bel] and [Nv] for the isotropic inverse problems. As we will explain below it is 
important to obtain nniqneness with mnch less regnlar coefficients. 

Once nniqneness resnlts have been obtained, one has to face np to the issne 
of ill posedness of the inverse problem. That is, one can make large changes in 
M that have only small effects on bonndary data obtained from examining the 
bonndary behavior of the resolvent kernel mentioned above. For example, given 
{M,g), one conld take an anxiliary manifold X, withont bonndary, of the same 
dimension as M, remove a small ball from X and from the interior of M, and 
connect these manifolds by a thin tnbe. One is faced with the task of stabilizing 
this ill posed inverse problem. One ingredient in this process involves having some 
a priori knowledge of the qnantities one is trying to determine, typically expressed 
in terms of a priori bonnds on these qnantities in certain norms. 

An early resnlt in this direction for the Gel’fand problem was given in [Al] by 
G. Alessandrini, who obtained conditional stability for the operator dive grad in a 
bonnded domain in Enclidean space, where e is a positive fnnction (scalar condnc- 
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tivity), assumed to be bounded in some Sobolev space with s > 0. See 

also [StU] for a related result for an anisotropic metric tensor close to Euclidean. 
Despite these successes, there is a clear need for coordinate-invariant constraints. 

In the case of trying to determine an unknown Riemannian manifold with bound¬ 
ary M, from boundary spectral data, it is natural to make a priori hypotheses on 
geometrical properties of M. Furthermore, if one must make such a priori hypothe¬ 
ses, it is desirable to get by with as weak a set of hypotheses as possible. There is 
then a tension between the desire to make weak a priori hypotheses and the need 
to establish uniqueness results. (For preliminary results in this direction see [K2L], 
[KaJ.) _ 

Here we impose a priori sup norm bounds on the Ricci tensor of M, and of dM. 
This, together with a Lipschitz norm bound on the mean curvature of dM ^ M, is 
shown in §2 to imply certain regularity, up to the boundary, of the metric tensor of 
M, when one is in “boundary harmonic coordinates” (dehned in §2). To be precise, 
we obtain regularity in the Zygmund space C^{M), a degree of regularity better 
than for any r <2 and just slightly worse than . This result has the following 
important advantage over a estimate. The Hamiltonian vector held associated 
with the metric tensor has components with a log-Lipschitz modulus of continuity. 
Hence, by Osgood’s theorem, it generates a uniquely dehned geodesic how, on the 
interior of M, and also for geodesics issuing transversally from dM. This property 
will be very important in §4. (We note that in the context of diherential geometry 
Zygmund-type spaces go back to the habilitation thesis of B. Riemann.) 

In §3 we obtain a compactness result for families of compact Riemannian mani¬ 
folds, of dimension n, with boundary, for which there are hxed bounds on the sup 
norms of RicM and Ric^M, on the Lipschitz norm of the mean curvature of dM, and 
on the diameter, and hxed lower bounds on the injectivity and boundary injectivity 
radius. We show that a sequence of such Riemannian manifolds has a subsequence, 
converging in the C'”-topology, for all r < 2, whose limit (M, ^f) has metric tensor 
\nCl{M). 

In §4 we study Gel’fand’s inverse boundary problem, recast in the form of an 
inverse boundary spectral problem. We show that, having boundary spectral data, 
we can recognize whether a given function h G C{dM) has the form h{z) = r^iz) = 
dist(a;, z), for some x G M, all z G dM, thus recovering the image in C{M) of M 
under the boundary distance representation. Such a representation, whose use was 
initiated in [Ku] and [KuL], plays an important role in the uniqueness proof, but 
for it to work we need to know that geodesics from points in dM, pointing normal 
to the boundary, are uniquely dehned. As noted above, this holds when the metric 
tensor in in C'^(M), and we obtain a uniqueness result in this category. This hts in 
perfectly with the compactness result of §3, to yield a result on stabilization of this 
inverse problem. 

Section 5 is devoted to the proof of several elliptic regularity results, of an ap¬ 
parently non-standard nature, needed for some of the hner results of §2. 
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Remark. A number of classes of function spaces arise naturally in our analysis. 
These include spaces C^{M), mentioned above. Here, if r = k + a, k E Z"*", a G 
(0,1), C'^{M) consists of functions whose derivatives of order k satisfy a Holder 
condition, with exponent a. The Zygmund spaces Cl{M) coincide with C^i^M) for 
r G (0, cx)) \ and form a complex interpolation scale. We also encounter L^- 
Sobolev spaces, and Besov spaces Bp p{dM), and bmo(M), the localized 

space of functions of bounded mean oscillation. Basic material on these spaces can 
be found in [Trl] and Chapters 2, 3 of [Tr2], in Chapter 13 of [Tl], and in Chapter 
1 of [T2]. 
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ported by NSF grant DMS-0072591. A. Katsuda’s work was partially supported by 
a Grant in Aid for Scientihc Research (C)(2) No. 14540081 of JSPS. M. Lassas’s 
work was partially funded by the Academy of Finland. M. Taylor’s work was par¬ 
tially supported by NSF grant DMS-0139726. Part of this work was done at MSRI 
and at the Oberwolfach RiP Program, and their support is gratefully acknowledged. 


2. Boundary Regularity for the Ricci Equation 

In this section we establish the key results on local regularity at the boundary of 
a metric tensor on which there are Ricci curvature bounds and a Lipschitz bound 
on the mean curvature. Our set-up is the following. 

Let R be a ball about 0 G M"", = B D {x : x'^ > 0}. Let T, = B (1 {x : x'^ = 0} 

and set O = U E. Let be a metric tensor on O, and denote by h its restriction 
to E. We make the following hypotheses: 

(2.0.1) Qjk G i7^’^(0), for some p > n , 

(2.0.2) G R^’2(E), l<j,/c<n-l, 

(2.0.3) Ric“ G L°°(0), 

(2.0.4) Ric^ G L°°(E), 

(2.0.5) H G Lip(E), 

Here H denotes the mean curvature of E C O, i.e., H = Tr A/{n — 1), where A 
is the Weingarten map, a section of End(TE). Our goal is to establish the following 
result. 

Theorem 2.1. Under the hypotheses (2.0.1)-(2.0.5), given z G E, there exist local 
harmonic coordinates on a neighborhood U of z in Ul with respect to which 

( 2 . 0 . 6 ) gjkeCliU). 
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Here C‘l{U) is a Zygmund space, as mentioned in §1. The harmonic coordinates 
for which (2.0.6) holds are arbitrary coordinates {v }^..., u'^) satisfying /S.u^ = 0 on 
a chart not intersecting E. On a neighborhood of a point in E, these coordinates 
are “boundary harmonic coordinates,” which are dehned as follows. We require 
(u^,...to be dehned and regular of class at least on a neighborhood of z 
in bl, and = 0. We require that = u^\t, be harmonic on E, i.e., annihilated 
by the Laplace-Beltrami operator of E with its induced metric tensor. We require 
to vanish on E, and we require (u^,... ,u^) to map a neighborhood of ^ in O 
diffeomorphically onto O. 

Let us note that the hypotheses (2.0.1)-(2.0.2) imply that various curvature 
tensors are well dehned. If {gjk) is the n x n matrix representation of a metric 
tensor in a coordinate system, its matrix inverse, the connection 1-form T is 

given by 

(2.0.7) T‘^bj = ^g'^'^idjgbm + dbgjm - dmgtj)- 

The Riemann tensor is then given by 

(2.0.8) 77 = dr + r A r. 

It is a matrix valued 2-form with components R°'bjk- We see that 

0 91 ^ ^ 77 ^’"( 0 ) ^ r G l2(0), R\jk e + L\n) 

^Ricbk e +L\n). 

The hypothesis (2.0.1) is stronger than the hypothesis in (2.0.9). It implies gjk G 
C'’’(0) for some r > 0, so (2.0.9) is applicable both to gjk on O and, in view of 
(2.0.2), to hjk on E. Furthermore, 

(2.0.10) gjk G 77i’P(0), p>n^ R\jk, Ric.fc, Ric^' fc, ^ G ^/-'’^(O), 

where Ric^ k = g^^ Ric^A: and S = Rie^ j is the scalar curvature of O. We mention 
parenthetically that one can use the fact that pointwise multiplication gives a map 

(2.0.11) X 

n — 1 

to obtain 

(2.0.12) gjk G (7(0) n 77^’2(0) ^ Ric^' k, S e 

However, we will not make use of (2.0.12) here. 
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We next consider the implication of (2.0.1) for the Weingarten map associated 
to E ^ O. The nnit normal to E is a vector held with coefficients 

(2.0.13) 

which by the trace theorem belongs to the Besov space i?p,p^^^(E). It follows that 
the Weingarten map has the property 

(2.0,14) A e B-y’’(E), 

as a consequence of (2,0.1), Thus we have a priori that H e Bp,y^{12), and the 
hypothesis (2.0.5) strengthens this condition on H, in a fashion that is natnral for 
the desired conclnsion of Theorem 2.1. 

Onr approach to the proof of Theorem 2.1 is to obtain the resnlt as a regnlarity re¬ 
sult for an elliptic boundary problem. We use the PDE (1.0.1) (the “Ricci equation”) 
for the components of the metric tensor, in boundary harmonic coordinates, and use 
Dirichlet boundary conditions on some components of gjk and Neumann boundary 
conditions on complementary components; see (2.1.8) and (2.1.16)-(2.1.17) for a 
more precise description. 

We will approach the proof of Theorem 2.1 in stages. In §2.1 we prove that 
the conclusion (2.0.6) holds when the hypothesis (2.0.1) is strengthened to gjk € 
C'^+®(n), for some s > 0. In §2.2 we replace (2.0.1) by the hypothesis that gjk € 
for some p > 2n. In §2.3 we prove the full strength version of Theorem 
2.1. These stages serve to isolate three rather different types of arguments, each of 
which is needed to prove Theorem 2.1, but which are perhaps more digestible when 
presented separately. Section 2.4 has some complementary results on the degree of 
regularity of the harmonic coordinates mentioned in Theorem 2.1. In section 2.5 
we demonstrate the non-branching of geodesics for metric tensors satisfying (2.0.6), 
including geodesics starting at a boundary point, in a direction transversal to the 
boundary. We also discuss examples of branching geodesics, for metric tensors only 
mildly less regular than those of (2.0.6), extending some examples of [Ha]. 

The version of Theorem 2.1 established in §2.1 is already useful for the results of 
§§3-4, and the reader particularly interested in §§3-4 could skip §§2.2-2.3, on hrst 
reading. However, the hypothesis (2.0.1) has a “natural” quality that we believe 
makes the additional effort required to work with it worthwhile. The arguments in 
§§2.2-2.3 require several elliptic regularity results that do not seem to be standard 
in the literature that we know, and their proofs are collected later, in §5. 


§2.1: First regularity result 

Here we prove that the conclusion of Theorem 2.1 holds when the hypotheses 
(2.0.1)-(2.0.2) are strengthened a bit. 
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Proposition 2.1.1. In the setting of Theorem 2.1, replace hypotheses (2.0.1)- 
(2.0.2) by 

(2.1.1) Qjk G for some s E (0,1), 

and retain hypotheses (2.0.3)-(2.0.5). Then the conclusion (2.0.6) holds. 

To begin our demonstration, let h denote the metric tensor induced on E, with 
respect to which (2.0.4) holds. By (2.1.1), hjk E C'^+'^(E), so there exist local 
harmonic coordinates ..., v'^~^ on a neighborhood (9 of z in E. Now we can hnd 
harmonic ..., on a neighborhood of z in O such that on O. Also 

we can hnd u'^, harmonic in Tl, with 'u"^|e = 0 and arrange that dx^^u^iz) ^ 0. We 
will have 

(2.1.2) e C'2+"([7). 

We have du^{z ),..., du^{z) linearly independent, so, after perhaps further shrinking 
U we have a harmonic coordinate chart on U, a set we relabel as O. As mentioned 
below the statement of Theorem 2.1, this is what we call a set of boundary harmonic 
coordinates. In these new coordinates, (2.1.1) and (2.0.3)~(2.0.5) are preserved. 
Now in harmonic coordinates the metric tensor satishes the elliptic PDE 

(2.1.3) 

where A acts componentwise on gim, as 

(2.1.4) Au = g~^/^dj{g^/^g^'"dku), g = det{gjk), 
and 

(2.1.5) Fim = Bimig,Vg)-2Ric^^. 

Here Bim is a quadratic form in V (7 with coefRcients that are rational functions of 
gjk- Thus, from (2.1.1) and (2.0.3) we have 

( 2 . 1 . 6 ) FemeL^in), 

and the coefRcients of A have the same degree of regularity as gjk in (2.1.1). 

Now, if j, /c < n — 1, then well known local regularity results on E following from 
(2.0.4) give 

(2.1.7) gjk \,, = hjkeH ^’ P { Jl ), Vp<cx), 

but in fact there is the following rehnement of (2.1.7), established in Proposition 
III.10.2 of [T2]: 

(2.1.8) 9jk\j:^ = hjk E l<j:k<n-l. 

Here 1)2’°° denotes the bmo-Sobolev space of functions whose derivatives of order 
< 2 belong to bmo, the localized space of functions of bounded mean oscillation. 
We establish the following (after perhaps shrinking to a smaller neighborhood of 
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Lemma 2.1.2. Under our working hypotheses we have, in the harmonic coordinate 
system {u ^,..., u'^), 

(2.1.9) gjkeCUU), l<j,k<n-l. 


Proof. First, extend Fim by 0 on \ O and solve Awim = Fim on a neighborhood 
of 0 in B, where we obtain A in the form (2.1.4) with g^’^ extended across E to 
gjk g C'^+®(i3). Local elliptic regnlarity resnlts imply 

( 2 . 1 . 10 ) 

It follows that W£m\T. £ C'*(S) and, via (2.1.8), 

(2.1.11) gjk-Wjk\j. = bjkeCl{T,), j,k<n-l, 

while 


(2.1.12) A{gjk-Wjk)=0 on O. 

Given onr assnmed regnlarity of the coefficients of A, standard Schander resnlts 
give 

(2.1.13) bjk e ^(E) gjk - Wjk e ^’^(0), l<r<2, 2<r<2 + s. 

Actnally, the case 1 < r < 2 is perhaps not so classical, bnt see [Mol], Theorem 7.3 
or [GT], Gorollaries 8.35-8.36. From here, an interpolation argnment gives 

(2.1.14) bjk e CliT) ^ gjk - Wjk e C^m. 

See [Tl], Ghapter 13, §8, particnlarly (8.37), for interpolation in this context. This 
establishes (2.1.9). 

To continne, following [An2], we switch over to PDF for g^'^. Parallel to (2.1.3), 
we have 

(2.1.15) Ag^^ = B^^{g, Vg) + 2(Ric^)^"^ = F^^, 

and (2.1.1) and (2.0.3) give G L°°(0). We take m = n and proceed to derive 
Nenmann-type bonndary conditions for the components g^'^, 1 < i < n. In fact, as 
we will show. 


(2.1.16) Ng^^ = -2(n - l)Hg^^, on E, 

and, for 1 < £ < n — 1, 


Af/" = -{n- l)Hg”' + 




1 1 


tk fA ^nn 


(2.1.17) 


2 


on E. 
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Here H is the mean curvature of E, which we assume satishes (2.0.5), and N is the 
unit normal held to E, pointing inside O. 

To compute (2.1.16)-(2.1.17), we use 

V/ni'n 1 

(2.1.18) = (Vu^ Vu-), iV = ^ 
and 

(2.1.19) = {Vn^u\ Vu^) + Vn^u^)- 

We also use the fact that is harmonic on O and u^\y^ = is harmonic on E (0 
if £ = n). 

Note that if {ej : 1 < j < n — 1} is an orthonormal frame on E and X a vector 
held on fl (say both having coefficients in (7^(0)) then 

n—1 

(2.1.20) div = 5^(VeW,e,) + {VnX,N). 

i=i 

In particular, for Xi = Vu^, we have divW^ = Au^ = 0, so the hrst term on the 
right side of (2.1.19) is equal to —times 

n—1 

( 2 . 1 . 21 ) 

1 = 1 

Let us set 

(2.1.22) = Wf + Xj, Wf = {Xe, N)N = ipN, Xj = Vv\ 

with Xj tangent to dVt and ip = Since Vu^, e^) = divVu^ = 

Av^ = 0, we have (2.1.21) equal to 



so the hrst term on the right side of (2.1.19) is equal to —(n — l)Hg^'^. The case 
£ — n gives (2.1.16), since the two summands in (2.1.19) are then the same. 



10 


To continue when £ n, we note that 

(2.1.24) (Vu^VnVu^) = (iV, Vx, V 11 -), 
with — Vu^. In fact, generally a 1-form rj satishes 

dr^{X,Y) = {Vxr),Y)-{VYr^,X), 

and applying this to ry = du^ and X = Xi, Y = Xn gives (2.1.24). Now 

Xeg^^ = 2{Vx,Vu^,Vu^), 

so (2.1.24) is equal to 

(2.1.25) I _ (Xe, = I _ g^’^dkg'^^, 

which gives (2.1.17). 

Having (2.1.15)-(2.1.17), we can establish further regularity of the functions g^'^. 
Lemma 2.1.3. In the harmonic coordinate system ... ,u'^), we have 

(2.1.26) g^^eCliTt), l<i<n. 

Proof. As in Lemma 2.1.2, we extend g^^ to g^^ G and extend the right 

side of (2.1.15) by 0 and produce a solution 

(2.1.27) e [)2’“(H) C Cl{B) 

to Aw^'^ = . Then g^'^ — satishes 

(2.1.28) A(/" - w^") = 0, on O, 
and we have 

(2.1.29) Nig^^-w^n\^eCliY), 

by (2.1.16) and (2.1.27), plus the regularity N G C'^"*'^(E). As in Lemma 2.1.2, we 
can apply an interpolation argument to Schauder-type estimates (see §5.3) and get 

(2.1.30) 

and hence (2.1.26) holds for i — n. Having this, we get from (2.1.17) that 

(2.1.31) N{g^--w^^)\^eCl{Y), 

which gives g^'^ — G (7^(0), so (2.1.26) holds for all i. 

The hnal step is to verify that Lemmas 2.1.2 and 2.1.3 yield regularity of gin- 
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Lemma 2.1.4. In the setting of Lemmas 2.1.2-2.1.3, 

(2.1.32) gne^gen^C^A^- 

Proof. Let g = det{gjk) and let Aim be the determinant of the (n — 1) x (n — 1) 
matrix formed by omitting colnmn I and row m from the matrix {gjk)- Then 

(2.1.33) g^'^ = tUllA.k. 

g 

By Lemma 2.1.2, Ann ^ C'*(b^)- Applying Lemma 2.1.3 to g'^'^ (which is > 0) we 
have 

(2.1.34) g = AnJg^^eCl{Tl). 

Then it follows that 

(2.1.35) Atn = Ani^{-lT^^gg^^eCl(Sl), l<l<n. 

Another way of pntting this is the following. Let 

(2.1.36) hjk = gjk, l<j,k<n-l; h = det{hjk) 

and (h^^) be the matrix inverse to (hjk). Then A^nA ^ n — 1, can be written in 
the form 

(2.1.37) Atn = {-ir-^+^gjnhh^^. 

Now the regnlarity and positive-dehniteness of {hjk)i<j,k<n-i applied to (2.1.37) 
yield 

(2.1.38) igni,...,gnn-i)eCl(^)- 
Finally, the identity 

(2.1.39) gjng^^ = 1, 

the regnlarity of g^'^ in (2.1.26) and of gjn for j < n — 1 in (2.1.38), pins the fact 
that g'^'^ > 0, yield 

(2.1.40) gnneCA^), 

proving the lemma, and completing the proof of Proposition 2.1.1. 
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§2.2: First improvement 

In this section and the next we obtain regnlarity with a weaker a priori hypothesis 
than (2.1.1). As we noted above, the results of §2.1 suffice for the applications in §3, 
but these improvements are quite natural (if not trivial to implement) and surely 
have the potential for applications elsewhere. 

Here we do strengthen the hypothesis (2.0.1) to some degree. Namely we assume: 

(2.2.1) G H'^’^(O), p > 2n, 

We retain hypothesis (2.0.2), i.e., 

(2.2.2) hjfc G H'^’^(S), l<j,k<n-l. 

Here Q and E are as in §1 and hjk = Qjkli:, for 1 < j,k < n — 1. Note that (2.2.1) 
implies gjk G (7^(0) with r = 1 — n/p > 0, and hence hjk G C'’'(E). 

Our strategy is to show that the hypotheses (2.2.1)-(2.2.2) together with (2.0.3)- 
(2.0.5) imply gjk G C'^"*'®(fl) for some s > 0, so Proposition 2.1.1 applies. In fact, 
we will expand the scope of the investigation here, and establish this conclusion 
under the following hypotheses, which are weaker than (2.0.3)-(2.0.5): 


(2.2.3) 

Ric^ G 

Pi > n, 

(2.2.4) 

Ric^ G LP^E), 

P 2 > n - 1, 

(2.2.5) 

H G (^"(E), 

s > 0. 


Our goal in this section is to prove: 

Proposition 2.2.1. Under the hypotheses (2.2.1)-(2.2.5), given z G E, there exist 
local harmonic coordinates on a neighborhood U of z in Q with respect to which 

( 2 . 2 . 6 ) gjkeC^+^{U), 

for some s > 0. 

As before, we begin by constructing local harmonic coordinates on a 

neighborhood O of z in E. Knowing that hjk G C^{E), we can do this, and making 
use also of hypothesis (2.2.2) we have 

(2.2.7) G C^+^O) n 

by Proposition 9.4 in Chapter HI of [T2]. It follows that (2.2.2) persists in this new 
coordinate system. As a consequence of the fact that hjk G i?p,p^^^(E), we also 
have 


( 2 . 2 . 8 ) 
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This result is established in §5.1. 

Next we hnd harmonic functions v }^on a neighborhood U oi z in Vt 
such that on O and we hnd harmonic in O, with 'u"^|e = 0 and arrange 

that dx^u^{z) 7 ^ 0. Given (2.2.7)-(2.2.8) and the hypothesis (2.2.1), we claim that 

( 2 . 2 . 9 ) ..., e C^+^(U) n 

The fact that G follows from Corollaries 8.35-8.36 of [GT], or [Mol], 

Theorem 7.3. The fact that G i7^’P(t7) is established in §5.2. We have du^{z ),..., 
du^{z) linearly independent, so after perhaps further shrinking U we have a har¬ 
monic coordinate chart on 77, which we relabel O. In these new coordinates, (2.2.1)- 
(2.2.5) are preserved. 

In fact, now that we have switched to harmonic coordinates, we can improve 
(2.2.2), making use of (2.2.4). It follows from Proposition 10.1 in Chapter III of 
[T2] that 

(2.2.10) hjk e l<j,k<n-l. 

In particular, 

(2.2.11) hjk G G^^'^(s), for some s > 0. 

We may as well suppose s G (0,r). Now we can prove: 

Lemma 2.2.2. In the harmonic coordinate system (u^,... ,u^), 

(2.2.12) 1 <i,m<n-l. 

Proof. We know girn solves the Dirichlet problem 

(2.2.13) ^gi^ui Firm dim \ hi.^., 1^-7, m F n 1, 

where 

(2.2.14) Fim = Bimig^Vg) - 2Ric2n ■ 

From (2.2.1) and (2.2.3) we have 

(2.2.15) Fgm e Qi = min(p/2,pi) > n, 

and the coefficients = g^^‘^g^^ of A are known to belong to C G^(f2). 

The next step in the proof is by now familiar. Extend Fi^ by 0 on \ O and 
solve Avim = Fim on a neighborhood V of ^ in O with Vim G G^+^CP), s > 0. Then 
Wim = gim - Vim SolveS 

Awim = 0 on O, Wim\j. = him “ Vim\j^ € G^+^(S), 


( 2 . 2 . 16 ) 
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and the previously cited results of [Mol] and [GT] yield wirn ^ hence 

( 2 . 2 . 12 ). 

It remains to show that 

(2.2.17) €^1+^(0), l<£<n, 
with s > 0. In fact, if we show that 

(2.2.18) e G^+"(n), 1 < £ < n, 

then an argument parallel to the proof of Lemma 2.1.4 yields (2.2.17). 

As before, we have 

(2.2.19) A/” = Vg) + 2(Ric^)^^ = 

As in (2.2.15), we have 

(2.2.20) G L'^^(O), qi = min(p/2,pi) > n. 

However, this time it is not so straightforward to produce the Neumann-type bound¬ 
ary conditions (2.1.16)-(2.1.17). 

Consider (2.1.16). The right side is well dehned; we have Hg'^'^\Y, G G'^(E), for 
some s > 0. As for A, the unit normal held to E is also Holder continuous of 
class . But applying N to g'^'^ G i7^’P(0) does not yield an object that can be 
evaluated on E. One has the same problem with the left side of (2.1.17), and the 
right side of (2.1.17) is also problematic. 

Our next goal is to show that a weak formulation of the Neumann boundary 
condition is applicable. Generally, the weak formulation of 

(2.2.21) Aw = F, Nw\^ = G 

is that for all test functions '0, i.e., all '0 G G°°(0) with compact support (inter¬ 
secting E but not the rest of 90), 


( 2 . 2 . 22 ) 


(Vm;, dV = - Fi^dV- dS. 


f2 


n 


Here dV is the volume element on O and dS the area element on E, both determined 
by the metric tensor in the usual fashion. Note that the left side of (2.2.22) is well 
dehned for all test functions '0 whenever Vrc G L^(0) and the right side of (2.2.22) 
is well dehned whenever F G L^(0) and G G L^(E). 
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Lemma 2.2.3. The function w = ( 7 "'” satisfies (2.2.22), with 

F = e qi>n, 

(2 2 23') 

^ G =-2(n-l)Hg^^l^eC^(T), r > 0. 

Hence the result that G for some p > 2n is improved to 

(2.2.24) g^^eC^+^m, s > 0. 


That (2.2.22) holds in this context seems a natnral generalization of (2.1.16), bnt 
the proof, given below, reqnires some work. Given that (2.2.22) holds, the regnlarity 
resnlt (2.2.24) follows from the results §5.3. 

Once Lemma 2.2.3 is established, we see that the right side of (2.1.17) is well 
dehned, and we can formulate: 

Lemma 2.2.4. For 1 < i < n — 1, the function w = g^'^ satisfies (2.2.22), with 


(2.2.25) 


F = G L'?i(0), qi>n, 

G = -{n- l)Ha>% + 6 C’-(E), 


r > 0. 


Hence the result g^'^ G i7^’P(0), for some p > 2n, is improved to 


(2.2.26) 


/"GGi+"(n), s>0. 


To set up the proof of Lemmas 2.2.3-2.3.4, let Oc = { 3 :^ £ ^ u^{x) > c} C O, 

for small c > 0, and let Ec = {x G O : u^(x) = c}. Since u'^ G the surfaces Ec 
are uniformly G^"'“®-smooth. Parallel to (2.2.10), we have from (2.2.3) that 

(2.2.27) 

so N is well dehned for small c > 0. Calculations parallel to (2.1.18)-(2.1.23) 
give 

(2.2.28) = -2(n - l)H,g^^, 

where H^ denotes the mean curvature of E^, i.e., (n — 1)77^ = Tr Ac, where 
denotes the Weingarten map of Eg. Note that, for X, Y tangent to Ec, 

(2.2.29) {A,X,Y) = {VxN,Y), N = {yu^y=gJ^. 


The assumption (2.2.1) implies 

(2.2.30) iV e A, e 
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For a fixed c > 0, from (2.2.27), we have N\y,^ G (Ec), hence Ac G 

(Ec), but (2.2.30) holds uniformly as c —0, and if E^ is identihed with 
E via the {u^, ■ ■ ■ , 'u"^“^)-coordinates, we have Ac continuous in c as c —0, in the 
space ^(E). 

Now pick a test function ip. From (2.2.28) we have, for each (small) c > 0, 


(2.2.31) 


j Vip) dV = - j F^^iP dV- j Gctp dS, 


where F'^'^^Gc are as in (2.2.23) with Ec instead of E. We let c —0. Since we 
already have 'Vg^'^ G L^(f2), the left side of (2.2.31) converges to 


(2.2.32) 


j{Vg^-,ViP)dV, 

12 


and the hrst term on the right side of (2.2.31) converges to 


(2.2.33) 


J F^^iPdV. 

Q 


Finally, from (2.2.30) we have that He —^ Hq in with p > 2n, which, 

given our knowledge at this point that gA g i7^’^(0), is more than enough to 
imply that the last term in (2.2.31) converges to 


(2.2.34) 


/ 


GiPdS. 


This proves Lemma 2.2.3. 

To proceed, we have, for 1 < £ < n — 1, c > 0, by a calculation parallel to 
(2.1.18)-(2.1.25), 


(2.2.35) 



- (n - l)i7c/” - 
1 1 

I ^ ^nn\ 


Here is the Laplace operator on the surface Ec, with its induced Riemannian 
metric tensor. Hence, given a test function ip^ we have, for each c > 0, 


f{Vg^-,ViP)dV 

fie 



GiP dS, 


fic 




(2.2.36) 
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with as in (2.2.25) and Gc given by the right side of (2.2.35). Again we let 
c —0, and since we know G L^(0), the left side of (2.2.36) converges to 

(2.2.37) j {Vg^^,V'il;)dV. 

ft 

Again the hrst term on the right side of (2.2.36) converges to 


(2.2.38) 


j F^^tljdV. 
n 


The last term in (2.2.36) is eqnal to 

(n-1) J Heg^^tl;dS + J 

/ \ Sf’ Sr 

(2.2.39) ^ 

- - / -^=g^’"dkg""tl^dS. 

2 7 

Sc 


As c —0, the hrst term in (2.2.39) converges to 


(2.2.40) 


(n 


1 ) j Hg^^i;dS, 

E 


by the same argnments as above (via (2.2.30)). Next (making nse of (2.2.1)) we 
have, for some p > 2n, 


(2.2.41) G 


u 




p,p 






with nniform bonnds and convergence as c —0, so the second term in (2.2.39) 
converges to 0 as c —0. 

Finally, since we already have g^^ G C'^"'“®(n), pins g^^ G the conver¬ 

gence of the last term in (2.2.39) as c —^ 0 follows, so Lemma 2.2.4 is proven. 


§2.3: Proof of Theorem 2.1 

In this section we hnish the proof of Theorem 2.1 assnming that n < p < 2n. 
Going fnrther, we extend Proposition 2.2.1 as follows. 
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Theorem 2.3.1. Replace hypothesis (2.2.1) by 

(2.3.1) Qjk € n < p < 2n, 

and retain hypotheses (2.2.2)-(2.2.5). Then, given z G E, there exist local harmonic 
coordinates on a neighborhood U of z in Q with respect to which 

(2.3.2) gjkeC^^^{U) 


for some s > 0. 

To begin, we have local harmonic coordinates v^,..., v'^~^ on a neighborhood O 
of z in E satisfying (2.2.7)-(2.2.8) and then local harmonic coordinates ... ,u'^ 
as in §2.2, satisfying (2.2.9), and in these new coordinates (2.3.1) and (2.2.2)-(2.2.5) 
are preserved. We also continne to have (2.2.10)~(2.2.11). We next establish the 
variant of Lemma 2.2.2 that holds in this context. 

Lemma 2.3.2. In the harmonic coordinate system (n.^,... ,'u"'), we have 


(2.3.3) 
with 

(2.3.4) 


G iL^’"'i(0), 1<£, m<n-l. 


ri 


Qi 

1 - qi/n 


j. n p 

for any - < qi < -■ 


Proof. We continne to have (2.2.13)-(2.2.15), except that in (2.2.15) no longer have 
qi > n (rather n/2 < qi < n). As for the regnlarity of gem\j:, we continne to have 
(2.2.11). Extend Firn by 0 on \ O and solve Av£rn = on a neighborhood 
V of z in B. This time we can say V£m G iL^’^i(V), with ri as in (2.3.4). Since 
L'3'i c this follows from Proposition 1.10 in [T2], Chapter III. Then W£rn — 

him — V£rn satishes (2.2.16), so as before we have W£m G iL^’'^i(0), and this gives 

(2.3.3). Note that 

Tl 

(2.3.5) -^ri> 2qi. 

To proceed, we note that the resnlts (2.2.27)-(2.2.34) hold nnder onr relaxed 
hypotheses on p, so (2.2.23) holds, except we have qi as in (2.3.4) instead of qi > n. 
With this, we can prove the following. 

Lemma 2.3.3. For all e > 0, we have 

(2.3.6) g^^ G 


Proof. We hrst extend F'^'^ by 0 on \ O and solve 

(2.3.7) = gO^ e 
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on a neighborhood V of 2 ; in B, as in the proof of Lemma 2.3.2. In fact, by Propo¬ 
sition 5.2.2, we have 

(2.3.8) G 

since l/p + 1/p < 1/qi in this case. Of conrse (2.3.8) implies (2.3.7), bnt it also 
implies 

(2.3.9) s, = 

Ti 

the latter inclnsion holding for all £ > 0. 

Hence = (/g” -|- gi'^ where 

(2.3.10) Ag^^ = 0, = G- Ngo"^. 

Here G = -2(n - as in (2.2.23), so we know 

(2.3.11) GgC'^(E), Ng//^ e 

for all £ > 0. It follows from Proposition 5.5.2 that 

(2.3.12) g^^ e V£>0, 

which gives (2.3.6). For nse below we also record the non-tangential maximal fnnc- 
tion estimate 

(2.3.13) (V^r)* e L"'“^(E), 

also established in §5.5. The meaning of the left side of (2.3.13) is the following. 
First, we have Vgi"^ continnons on the interior of fl. Next, for x E T,, 

(Vgrnx) = snp \Vgr{y)\, 

where F^, = {y G 11 : d{y, x) < 2 d{y, E)}. 

Lemma 2.3.4. For 1 < i < n — 1, and for all e > 0, we have 

(2.3.14) G 

Proof. Using (2.3.8) and (2.3.13), we can extend the analysis in (2.2.35)-(2.2.41), 
to conclnde that g^'^ is a weak solntion to 

(2.3.15) Ag^^ = Ng^^ = G, 
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with G as in (2.2.25), i.e., 

(2,3.16) G=-(n- m'/X + 


From what we know so far, we have 

(2.3.17) G'eL"i-^(E), V£>0. 

As in the proof of Lemma 2.3.3, we can then write 

(2.3.18) g‘'' = C + X, 
with 

(2.3.19) = 0, Ng[^ = G- e L"i-®(E), 
which by an analysis parallel to that of (2.3.10)-(2.3.11) gives 

(2.3.20) g[^ e 
and proves (2.3.14). 

Now an argnment parallel to the proof of Lemma 2.1.4 gives 

(2.3.21) e 77^’"i-"(0), V£>0, 
for all j, k < n. Now, for £ small and qi close to p/2, 

(2.3.22) n — £ > p, 

an improvement over the hypothesis (2.3.1), as long as p < 2n. Thns replacing 
(2.3.1) by (2.3.22) and iterating this argnment a hnite nnmber of times, we establish 
that actnally (2.2.1) holds. 

This proves Theorem 2.3.1. It also rednces Theorem 2.1 to Proposition 2.1.1, 
and hence proves Theorem 2.1. 

Remark. Theorem 2.1 remains valid if we change conditions (2.0.3)-(2.0.5) into 

(2.3.23) Ric^ebmo(O), Ric^ G bmo(E), HeGl{E). 
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§2.4: Complements on coordinates 

Under the hypotheses of Theorem 2.1, we know there are local bonndary har¬ 
monic coordinates with respect to which gjk G (7^(0). Here we show that, with 
respect to snch coordinates, any other bonndary harmonic coordinates are smooth 
of class so O has the strnctnre of a C^-manifold. Since the interior behavior 

is simpler to establish than the behavior at the bonndary, we conhne onr analysis 
to the following resnlt. 

Proposition 2.4.1. Assume the metric tensor on O is of class (7^(0). Let u G 
solve 

(2.4.1) Au = fe Cim, = h G CUdft). 

Then u G (7^(0) . 

Proof. Classical Schander estimates readily give u G for all s < 1, as in 

(2.1.2) ; we jnst need to go a little farther. The key to success is to abandon our 
former practice (which worked so well) of writing the Laplace-Beltrami operator A 
in divergence form, and instead write it in non-divergence form: 

(2.4.2) g^^djdkU = /i, u\q^ = h, 
where 

The hypothesis on gjk, plus the current handle we have on u, gives g~^^‘^dj{g^^'^g^’^) 
(dku) G Cl{Ll), hence fi G C^(0), under the hypothesis on / in (2.4.1). Now 
the hypothesis on g^^ is strong enough for standard Schander estimates to apply, 
yielding 

(2.4.3) forr G (0,1)U(1,2). 
Then an interpolation argument gives the conclusion stated in Proposition 2.4.1. 


Remark. Under the hypotheses of Theorem 2.1, we also have 

(2.4.4) g,keH^^P{n), Vp<oo, 

in boundary harmonic coordinates. Given this, we have, in parallel with (2.4.1), 

(2.4.5) Au = fe u\g^ = he ^ue 

for q G (1, cxd). The proofs are simple variants of arguments presented above. 
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2.5: Non-branching (and branching) of geodesics 


Suppose {gjk) is a metric tensor on an open set V C The geodesic equation 
can be written in Hamiltonian form as 

(2.5.1) 
where 

(2.5.2) G(xX)^\g‘Hx)iA. 

In other words, the geodesic flow is the flow on T*V = V x generated by 


(2.5.3) 


^ ^dG d dG d 


Osgood’s theorem states that a vector held generates a uniquely deflned flow pro¬ 
vided its coefficients have a modulus of continuity u){t) satisfying 


(2.5.4) 


•‘-'2 dt 


uj(t) 


= oo. 


See, e.g.. Chapter 1 of [Tl] for a proof. An example for which (2.5.4) holds is 


(2.5.5) 


u{t) =tlogj, 


which is just a bit rougher than the Lipschitz modulus of continuity. In fact, one 
has 


(2.5.6) / e Cl(U) lf{x + y)- f{x)\< C||/||c; wdi/l), 

with u(t) given by (2.5.5). See, e.g., [T2], Chapter I, for a proof of this classi¬ 
cal result. This applies to (2.5.1)-(2.5.3) provided gjk € (7^(12), so we have the 
following. 

Proposition 2.5.1. Let (gjk) ^ C'*(V) be a metric tensor on an open set V C R"". 
Then the geodesic flow is locally uniquely defined. 

Corollary 2.5.2. With O as in Theorem 2.1, suppose gjk G (7^(0). Then the 
geodesic flow is locally uniquely defined when applied to any initial point [z, ^) with 
z E T, and v = v{^) G TzG pointing inside Vt, transversal to E. 

Here ^ and v are related by = g^^{z)^k- To prove the corollary, let V be a 
collar neighborhood of O. One can extend gjk to gjk G (7^(12), and we have that 
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the geodesic flow is locally uniquely deflned on T*V. If the initial point is ( 2 , 0? 
described above, the geodesic has tangent vector v at 2 ;, so under our hypotheses 
the geodesic initially moves into O, for small positive t. This behavior is hence 
independent of the chosen extension of gjk to V. 

We now discuss some examples of metric tensors only slightly rougher than 
treated in Proposition 2.5.1, for which there is branching of geodesics. These ex¬ 
amples are variants of some produced by P. Hartman in [Ha]. As in [Ha], we take 

(2.5.7) ds^ = h{v){dv? + dv'^). 


As noted there, curves of the form v — v{u) are (variable speed) geodesics provided 


(2.5.8) 2^ = [ 1 +^H'{v), H{v) = \ogh{v). 

Multiplying by dv/du yields 


(2.5.9) 
hence 

(2.5.10) 


d i'dv\‘^ 
du\duJ 




dv 

du 




Thus V = v{u) solves (2.5.8) provided dv/du 7 ^ 0 and 


(2.5.11) />(«(«» = l+([|)f 

If, however, 

(2.5.12) h{0) = 1, A(0) = 0, 


then n = 0 solves (2.5.8), so if we produce another solution to (2.5.8) such that 
n( 0 ) = 0 , n'( 0 ) = 0 , we will have branching of geodesics. 

For the first class of examples, we pick k E {1, 2, 3,...} and construct h{v) such 
that 


(2.5.13) v{u)^u^^+^ 
solves (2.5.8). Indeed, by (2.5.11), this happens when 

(2.5.14) h{v) = l + i2k + i)2|^|4fc/(2fc+i)_ 


In this case the metric tensor has one derivative in k = 1 

was explicitly mentioned in [Ha]. 



24 


For another family of examples, we take k G {1, 2,3,...} and constrnct h{v) snch 
that 


(2.5.15) 


v{u) — e 


sgnn, 


solves (2.5.8). This time we have 
(2.5.16) 

and hence (2.5.11) gives 


u = 


sgn V fdv'^ 

2 p 2 

1 - ^ p-s/Gl" 

|log 

n \du/ 

' . y ^ 2 fc -|-2 


(2.5.17) h{v) = 1 + k'^l log 

For these examples, the hrst derivative of the metric tensor has a modnlns of continn- 
ity only slightly worse than log-Lipschitz. Also the metric tensor has two derivatives 
in for all p < cx). Hence the curvatnre belongs to for all p < cx). 


3. Geometric convergence for manifolds with boundary 

A sequence {Mk,gk) of compact Riemannian manifolds with boundary dM^ is 
said to converge in the C'^'-topology (given 0 < r < cx) to a compact Riemannian 
manifold (M, ^f) provided that is a C'’’ metric tensor on M and, for k sufficiently 
large, there exist diffeomorphisms : M ^ Mk such that F^gk converges to g in 
the C'^'-topology. (Necessarily Fk ■ dM dMk.) In this section we will identify 
classes of Riemannian manifolds with boundary that are pre-compact in the C^- 
topology, for any given r < 2. 

We work with families of Riemannian manifolds with boundary of the following 
sort. Fix the dimension, n. Given RoAo,‘S'o, <^o ^ (0, x), denote by Al(RoAo, 5'o, do) 
the class of compact, connected, n-dimensional Riemannian manifolds with bound¬ 
ary {M,g), with smooth metric tensor, with the following four properties: 

(3.0.1) II RiCM ||l°°(M) < Roi II hiCOM \\L°°{dM) < -Ro, 

where Ric denotes the Ricci tensor. 

(3.0.2) iM > *0, idM > *0, ib > 2fo- 

Here Im denotes the injectivity radius of M, Iqm that of dM, and A the boundary 
injectivity radius of M. 


( 3 . 0 . 3 ) 


-R||Lip(9M) < So, 
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where H is the mean curvature of dM in M. 

(3.0.4) diam(M, 5 f) < do- 

We recall the concept of boundary injectivity radius, ib- It is the optimal quantity 
with the following property. Namely, there is a collar neighborhood C of dM in M 
and a (unique) function / G C'^(C) such that = 0, | V/| = 1, f{C) D [0, ib). 

With this, local coordinates (n^,..., on an open set in dM can be continued 

inside, as constant on the integral curves of V/, to produce, along with = /, a 
set of “boundary normal coordinates.” 

To further clarify the hrst part of (3.0.2), we mean that 

(3.0.5) Expp : Bp{0) M, 

where Bp{0) = {n G TpM : g{v,v) < p^}, is a diffeomorphism for p = Iq if 
dist(p, dM) > ig and it is a diffeomorphism for p = dist(p, dM) if dist(p, dM) < ig. 
The main goal in this section is to prove the following. 

Theorem 3.1. Given Rg, ig^ 5'o, dg G (0, cx)), M(Rg, ig, Sg, dg) is precompact in the 
-topology for each r < 2. In particular, any sequence {Mk, g^) in M.{Rg,ig, Sg,dg) 
has a subsequence that converges in the -topology to a limit (M, g). Furthermore, 
the metric tensor g belongs to C‘l{M). 

Such a result was established in [Anl] in the category of compact manifolds 
without boundary; subsequently there have been expositions in [HH] and in [Pe]. 
Our proof of Theorem 3.1 follows the structure of the argument in [Anl], with 
necessary modihcations to treat the case of nonempty boundary. In this regard 
the boundary regularity results of §2.1 play a major role. The part of the 
conclusion is also more precise than that noted in earlier results. This precision will 
be of major value in the application of Theorem 3.1 to results on inverse boundary 
spectral problems presented in §4. 

The proof of Theorem 3.1 involves a blow-up argument that takes us outside 
the category of compact manifolds, and it is useful to have the following notion of 
pointed convergence of a sequence (M^, gk^Pk), with distinguished points pk G M^. 
We say (Mjt, gk,Pk) converges to (M, g,p) in the pointed (P^-topology provided the 
following holds. For large k there exist pk < o^k, Pk Z' cxo, and compact Ofc C Mk 
and Vk C M, such that 

BpkiPk) C flfc C Bfj^^fp}-), Bp^{p) C Vfc C B(j^[p), 
and diffeomorphisms 

Fk '.Vk ^ llfc, Ajt : Vfc n dM —blfc n dMk, 

such that F^gk converges to g in the C'^'-topology (on each compact subset of M) 
and F^^{pk) —^ p. We assume M is connected, and these hypotheses imply {M,g) 
must be complete. 
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§3.1: Basic convergence results 

We begin by describing a slight variant of a well known “abstract” convergence 
resnlt. As before, fix the dimension n. Given s = i + a (i G 0 < cr < 1), 
p G (0, cx)), Q G (1,2), let Af{s,p,Q) denote the class of connected, n-dimensional 
Riemannian manifolds (M, g) with bonndary, with the following properties. Take 
p G M. 

(i) If dist (p, dM) > p, there is a neighborhood t/ of p in and a coordinate 

chart (p : i?p/ 2 ( 0 ) —G, snch that, in these coordinates 

(3.1.1) Q~^\vf < 9jk{x)r]^r]’" < Q^\vf, 
and 

(3.1.2) p^ '^sup\x-y\~^\df^gjk{x)-d^gjk{y)\<Q-l. 

\p\=i 


(ii) If dist (p, dM) < p, there is a neighborhood of p in M and a coordinate chart 
(f : B^p{0) U snch that {x'^ = 0} maps to dM and (3.1.1)~(3.1.2) hold in these 
coordinates. 

Let p, Q) denote the class of pointed manifolds {M,g,p), with p G M, satis¬ 

fying these properties. 

The following compactness resnlt goes back to [Ch]; a detailed proof is given on 
pp. 293-296 of [Pe] for the case of manifolds withont bonndary, bnt no essential 
changes are reqnired for the case of manifolds with bonndary. 

Theorem 3.1.1. Given s, p E (0, cx)), Q G (1,2), the class p, Q) is compact 
in A/”* (s', p, Q) in the pointed -topology, for all s' < s. 

As for convergence withont specifying base points, we can dehne J\f{s, p, Q, do) 
to consist of (M, g) G J\f{s, p, Q) satisfying also diam (M, g) < do, and conclnde 
that Af{s, p, Q, do) is compact in A/’(s', p, Q, do) for ah s' < s. 

To apply Theorem 3.1.1 to onr sitnation, we will show that, given Ro,io, So,do G 
(0, cx), and given Q G (1, 2), s G (1, 2), there exists p > 0 snch that 

(3.1.3) M{Ro,io,So,do) C N'{s,p,Q). 

In fact, we will establish a result that is more precise, in two respects. 

For one, we will produce harmonic coordinates, in case (i), and boundary har¬ 
monic coordinates, in case (ii). Doing this brings in the notion of G^-harmonic 
radius, introduced in [Anl] in the context of manifolds without boundary. To be 
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precise, if s = £ + a, as above, a number = r^(p, g, Q) is called the C'^-harmonic 
radius of (M, g) ai p E M provided it is the optimal quantity with the following 
property. Take any p < rf^. Then, in case (i), there exist harmonic coordinates 
: U i?p/2(0) such that (3.1.1)-(3.1.2) hold on i?p/2(0). In case (ii), there 
exist harmonic coordinates : U such that (3.1.1)-(3.1.2) hold on 

5^(0). If M is compact, we dehne the C'^-harmonic radius of {M,g) by 

(3.1.4) rl{M,g,Q)= in^rl{p,g,Q). 

peM 

The containment (3.1.3) follows from the fact that, given Ro^io, Sq, do E (0, cx)) 
and s E (1,2), Q E (1,2), there is a lower bound on the C'^-harmonic radius of 
g) E M.{Ro: io: Sq, do). 

The result we will establish is more precise in one more respect; we bring in 
the notion of C'^-harmonic radius. Namely, a number th = rh{p, g, Q) is called the 
C'^-harmonic radius of (M, g) ai p E M given the circumstances described above 
(in the dehnition of C''^-harmonic radius) but with (3.1.2) replaced by 

(3.1.5) ^ s\\p\x-y\~^\d'^gjk{x)+d'^gjk{y)-‘id^gjk{{x + y)/2)\<Q-l. 
I/3|=1 


Then the C'^-harmonic radius of (M, g) is dehned by 

(3.1.6) rh{M,g,Q) = in^rh{p,g,Q). 

peM 

We note the following elementary but important scaling property of harmonic 
radius, valid for A G (0, cx)): 

(3.1.7) rh(M, X^g, Q) = Xrh(M. g, Q), r^(M, X^g, Q) = Xrl(M, g, Q). 

Our next goal is to show that, given Ro, io, So, do E (0, cx), there is a lower bound 
on the C'*-harmonic radius of (M, g) E M(Ro, io, So, do). 


§3.2: Harmonic radius estimate 

As advertised, our goal here is to prove the following. 

Theorem 3.2.1. Let Ro,io,So, and do he given, in (0, x), and let Q E (1,2) he 
given. Then there exists tm = xm^Ro, io. So, do, Q) > 0 such that 


(3.2,1) 


rh{M,g,Q) > tm, V {M,g) € io,SoVo). 
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The proof will be by contradiction. Snppose there exist (M^, gu) G ‘S'o, do) 

snch that 

(3.2.2) rh{Mk,hiQ) ^ ^ 0. 

Let ns scale the metric to gk = ^k'^gk: and consider the scaled Riemannian 
manifolds {Mk,gk)- Then 

(3.2.3) rhiMk,gk,Q) = I, 
while, for the rescaled metric, 

(3.2.4) II RiCfc ||l°°(M) 0, || RicgM*, ||L°°(9Mfe) 0, ||-f^A:||Lipi(9Mfc) Oj 

and 

(3.2.5) ik oo, idMk oo, ib,M^: oo, 

as /c —cx). We will show that conditions (3.2.3)-(3.2.5) lead to a contradiction. 

To see this, pick pk G Mk snch that 

(3.2.6) rh{pk,gk,Q) = I- 

By Theorem 3.1.1, there is a snbseqnence, which for simplicity we will also denote 
{Mk, gk,Pk) which converges to a pointed manifold {M,g,p) in the C'^'-topology, 
where r G (1,2) is arbitrary. With Tk = (iistk{x,dMk), where distfc is the distance 
on {Mk,gk), there are two possibilities: 

(i) Tfc(pfc) ^ CX). 

Then we claim {M,g) is isometric to with its standard flat metric. 

(ii) Tfc(pfc) <K<oc^ ^ 

Then we claim {M,g) is isometric to M_)_, with its standard flat metric. 

We demonstrate these claims in Lemma 3.2.2 below. For now we assnme this, 
and proceed with a farther analysis of these two cases. 

Case (i). This case is treated in [Anl]. We recall the argnment here, in a slightly 
varied form (also borrowing from [HH]), both to establish a slightly stronger con- 
clnsion and to set the stage to examine Case (ii). 

We have neighborhoods Uk of pk in Mk, identihed with = {x E M” : |a:| < 5}, 
with Pk = 0, and the metric tensors gk ^ S in C’^-norm on B^, where 6 is the 
standard Enclidean metric tensor on Taking {1 < u < n) to be the standard 
Cartesian coordinates on M"", we solve 

u’^ — x^ on dB^, 


(3.2.7) 


Akul = 0 on Rs, 
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where is the Laplace operator with respect to the metric tensor g^: 

(3.2.8) AkU = {9ll)9li)dju ). 

Note that 

(3.2.9) Afc('Ufc - x'') = on S 5 , ul - x’' = 0 on dB^, 
where 

(3-2.10) n = 

for which we have 

(3.2.11) IIAl’llc.-.(s.)^0, 

In view of this and the nniform estimates on the coefficients in the elliptic 
differential operators (3.2.8), elliptic regnlarity gives 

(3.2.12) ul-x’^ ^0 in k ^ 00 . 

Hence, for all large /c, (n.^,..., form a coordinate system on Uk, and if now 
denotes the components of the metric tensor gk in this coordinate system, then 

(3.2.13) Iblf - 52 ,||c^(B 4 ) ^0, /c ^ 00 . 

The last step is to obtain an analogne of (3.2.13) in a stronger norm, nsing the 
Ricci eqnation, which implies 

(3.2.14) At(4>-3.,)=f;y, 
where is as in (3.2.8), acting componentwise, and 

(3.2.15) = Bij{gira: Vgim) - 2 Ric(^^. 

Here Bij is a qnadratic form in ’Vgem- In view of (3.2.4) and (3.2.13) we have 
\\^iP\\L’=°(B 4 ) 0. Hence local elliptic regnlarity resnlts on (3.2.14) (pins another 

appeal to (3.2.13)) give 

(3.2.16) \\g^j^ - Sij\\c2jBs) ^0, /c ^ 00 , 
and even the stronger resnlt 


(3.2.17) 


00 . 
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Hence we have rhipk, 9kj Q) > 3 for large /c, contradicting (3.2.6). 

Remark. Clearly there is nothing special abont “3” in (3.2.16). For any L G (1, cxd) 
one obtains (3.2.16) with C^{Bs) replaced by C^{Bl). 

Case (ii). In this case, for large /c, there is a nniqne qk G dM^ closest to pk- 
In addition to the pointed convergence {Mk: Qk^Pk) —^ we also have 

{^Mk,gk,qk) {dM,g,q), and distkipk.qk) ^ dist(p, q ). This time {M,dM,g) = 
(R^, R"'“^, 5), with R"-"^ = {x"" = 0}, and we can identify q with 0 G R”"^ C R”. 
Fix L > 2dist(p, q) + 4. 

In this case we have neighborhoods Uk of qk in identihed with 

-Sl+5 = {x E R^ : |a:| < L + 5}, 

with qk = 0, and the metric tensors gk ^ S in C^-norm on Note that 

{x G Mk : distA;(a;,Pit) < 2} C B^^^^ for large k. 

With x'^ {1 < p <n) the standard Cartesian coordinates on R_|_ [x'^ = 0 dehning 
9R” ) we hrst solve, for 1 < z/ < n — 1, 

(3.2.18) ^dMkXk = 0 on Bl+ 5 , = x'^ on dBL+ 5 , 

where 

i?L+5 = {x G R"'”^ : |a:| < L + 5}. 

Parallel to (3.2.12), we have 

(3.2.19) k^O, l<u<n-l. 

Next solve 

(3.2.20) = 0 on R+^g, u'^k\d+Bi^+^ ^ ^ 

where 

9+Sl+ 5 = a(st+5) \ -Bl+5. 

Parallel to (3.2.9)-(3.2.11), we have 

(3.2.21) Z^kiK - x'^) = on R+^g, ul-x'^^ifl on ^R^+g, 
with 


(3.2.22) 


\\fk\\c-HB+^^) 


0 , 
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and 

(3.2.23) IIV^fcllL°°(9B+_l_g) 

This gives first a global estimate 

(3.2.24) — X ^ 0, 

and then a “local” estimate (away from the corners of 

(3.2.24) — X IIc’'+i(b+_|_ 4 ) 

as /c —> cx), for 1 < z/ < n — 1 . 

To constrnct we solve 

(3.2.26) /^kul = 0 on S+^ 5 , ul = rr" on 
Parallel to (3.2.25) we have 

(3.2.27) \\ul - ^0, ^ cx). 

Hence, for large /c, {u \^...,«.)!) form a coordinate system on with = 0} 

defining the face {x"^ — 0}. Parallel to (3.2.13), if denotes the components of 
the metric tensor in this coordinate system, we have 

(3.2.28) Wdij^ ~ ^b'llc’'( 9 B+_|_ 3 ) O’ /c —> 00 , 
for any r < 2 , 1 < t, j < n. 

As in Case (i), the last step is to obtain an analogne of (3.2.28) in a stronger 
norm. First, if we set = g^j^n ? 1 < b J ^ ^ ~ 1’ then the argnment in Case 

J J I _|_5 

(i) applies to yield 

(3-2'29) 

Now estimates snch as derived in §2 apply to glj \ solving (3.2.14)-(3.2.15) and the 
bonndary condition (3.2.29), for 1 < t, J < — 1, to yield 

(3.2.30) Wdij'^ ~ “^0, /c —> 00 , 

for 1 < t, J < — 1- Next, we nse 

(3.2.31) At(95, - «'") = Vs,^,) + 2 Ric[j), 
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with boundary conditions 

(3.2.32) ]V(9',",-3'’") = -(n-l)Ht 9 f^, + 1^^9f'|5j9Pl!;, 1 < t < n - 1, 

y^ik) 

on Sl+ 5 - The results of §2 together with hypotheses in (3.2.4) and the analogue of 
(3.2.28) for IIc-(b+_^ 3 )’ hrst 

(3.2.33) llfi'(fc) ~ ^ k ^ oo, 

and then 

(3.2.34) llfi'(A:) ~ O’ /c —>■ oo, i < n. 

Finally an argument parallel to that in Lemma 2.1.3 gives (3.2.30) for all i,j < n. 

Hence we have rh{pk, Qk, Q) > 2 for large /c, contradicting (3.2.6). This hnishes 
the proof of Theorem 3.2.1, modulo the proof of the following lemma. 

Lemma 3.2.2. Assume that a sequence of Riemannian manifolds {Mk, gk,Pk) sat¬ 
isfies the conditions of Theorem 3.2.1 and g^^Pk) —^ {M,g,p) in the C^- 
topology, r G (1,2). Assume that 

(3.2.35) II Ricfc llLoo(Mfc) —^ 0, || RicaM^ \\L°°{dMk) 0’ \\Hk\\L°°{dMk) 0’ 

and 

(3.2.36) iMk oo, idMk oo, ib,Mu oo, 

as k ^ oo. Then (M, g) is isometric to either M"' or M_|_. 

Proof. We consider separately the cases Tk{pk) —^ oo and Tk{pk) < K < oo. 

Case (i). Suppose Tk{pk) —^ oo. 

This case is effectively treated in [Anl]. We recall briefly the argument, since it 
also plays a role in Case (ii). In this case the hypotheses imply 

(3.2.37) RicM = 0, 

weakly. Hence the metric tensor g is smooth in local harmonic coordinates. Also 
any unit speed geodesic 7 (t) such that 7 ( 0 ) = p is dehned in M for all f G M. Take 
T G (0, cx)). We claim 7 is the shortest path from 7 (—T) to 7 (T). 

Consider B^Tip) = {x G M : dist(a:,p) < 6 T}. The hypotheses imply that there 
exists /co such that for all k > ko there are open sets in identihed with Bqt(j>) 
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via diffeomorphisms, such that Pk ^ P and (/a; —> (7 in {Bqt{p))- Also, for each 
X e B2t{p), i{x,gk) > 3T. 

Since (M, g) is smooth, there is some cq > 0 such that i(p, g) > cq. Say 7 ( 00 ) = 
q. We have unit speed geodesics 7 a; on {Mk,gk), dehned for \t\ < 2T, such that 
7A:(0) = p, 'Jkick) = q, Ck ^ Cq. By Arzela’s theorem there is a subsequence 
'jk cr, uniformly on t G [—2T, 2T]. We see that a is a geodesic on (M, g), a(0) = 
p, a(co) = q. Hence a = 'y. (It follows that the entire sequence 'jk —^ 7, not just a 
subsequence.) Since i{x,gk) > 3T for x G i? 2 T(p), we can deduce that 

(3.2.38) distfc( 7 fc(-T), 7 fc(T)) = 2T, 

for k > ko- Since "ik ^ 1 uniformly on [—T, T] and (/fc —(7 in (7°, the left side of 

(3.2.38) converges to dist( 7 (—T), 7 (T)). 

Thus (M, g) is complete and Ricci flat and each geodesic through p is globally 
length minimizing. It follows from the Cheeger-Gromoll splitting theorem [CG] that 
(M, g) is isometric to standard flat M"'. 

Gase (ii). Assume Tk{pk) < K < 00 . 

Let qk G dMk be the nearest point on dMk to pk- By (3.2.36) qk is uniquely 
dehned (for large k). Also {dMk, gk,qk) converges in the topology to {dM, g, q), 
and dist(p, q) = \im Tk{pk)- In this case the hypotheses yield 

(3.2.39) RicM = 0, RIc^m = 0, = 0. 

Hence g is smooth on M in local harmonic coordinates. 

We will make strong use of two equations. One is: 

(3.2.40) AfcTfc = Hk on {xk = c}. 

Here, in slight contrast to notation in §2, we set Hk equal to the trace of the 
Weingarten map, i.e., to n — 1 times the mean curvature of the surface {xk = c}. 
The other is (in boundary normal coordinates {z,Xk), z G dM): 

(3.2.41) drHk ^ -TrAl - Ric(k){dr, dr), 

where Ak is the Weingarten map of the surface {xk = c}. See [Pe], §§2.3-2.4. 
Under our hypotheses, there exist Sk ^ 0 such that 

(3.2.42) \\Ric{k){dT,dr)\\L<^ < Sk, \\Hk{0)\\L-- < Sk- 
Hence 

(3.2.43) drHk< - ^—Hl+Sk. 

n — 1 
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In particular, Hk{T, z) < Hkij) for 0 < r < ib,Mk: where 

(3.2.44) drUk ==-pyUfc + £fc, /^fc(0) = Sk- 

Hence 

(3.2.45) Hk{T, z) < niax(A/(n - l)ek,ek), 

since (3.2.44) forces //^(t) < 0 before Hkij) can be larger than the right side of 

(3.2.45) . 

We next estimate Hk (r, z) from below. Pick h > 0 and suppose 


(3.2.46) Hk{To,z) < —S, for some tq < z G dM. 

For short, we set Hkij) = Hki^r^z). If k is large enough that > 2(n - l)£fc, it 

follows that Hkij) < —5 for all r G [To,f6,Mfe), and hence 

(3.2.47) ^ r>T„. 

2[n — 1) 


As Hk{To) < —5, this implies 


(3.2.48) 


Hkir) < - 


2(n-l) 

To) - 2(n - l)/5' 


But this implies blow-up of Hkij) somewhere on r G [tq, To + 2(n — 1)/5], contradict¬ 
ing the fact that Hk^r) is hnite on r G [0, ib,Mk)i which contains [tq, tq + 2(n — 1)/S\ 
when 5 > 4(n — l)/ib,Mu- This contradiction shows that, for any given h > 0, 

(3.2.46) must fail for all sufficiently large k. Hence we have Sk ^ 0 such that 

(3.2.49) l-kfk(T)l<Sk for 0 < r < 

If t(x) = dist(T, dM), then t/j — > r in due to gk ^ ff in Using again that 
gk ^ g in together with (3.2.40), (3.2.49), the boundary condition Tk\dMu — 0, 
and Propositions 5.1.1 and 5.2.1, we deduce that, for all s < 2, n/(n — 1) < p < oo, 


(3.2.50) 


Tk^r in 


where pk ^ oo for k ^ oo. In particular, ib^M = +cxd. Applying again (3.2.40), 
(3.2.49) together with gk ^ g in (7^, we obtain that 

71 

(3.2.51) Hk^H in LP{B+{pk)), -r < P < oo, 

n — 1 
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i.e., 

(3.2.52) H{t) = 0, V T > 0. 

Now, parallel to (3.2.41), we have 

(3.2.53) drH -\-Ty = — KiCMidr^dr) = 

where A is the Weingarten map on the surface {r = c}, and hence Tr = 0, which 
implies 

(3.2.54) A = 0. 

Furthermore, since 

(3.2.55) drgrj =2A^gej =0, 
we have 

(3.2.56) gij{T, z) = gij{0, z), z e dM. 

By the argument of Case (i), dM is isometric to so this implies that M is 

isometric to [0, cx)) x dM = R_(_, with its standard flat metric. 

This finishes the proof of Lemma 3.2.2, hence of Theorem 3.2.1. 

§3.3: Proof of Theorem 3.1 

All the work needed to prove Theorem 3.1 has been done, and we need only collect 
the pieces. Say (M^, gt) G Ai(-Roj * 0 ? 'S'o, do)- The results of §3.2 show that Theorem 
3.1.1 is applicable. Hence, after passing to a subsequence, we have diffeomorphisms 
Ffc : M —Mfc and a Riemannian metric on M such that -F^gk g in C^{M) for 
all r < 2. It remains to show that g G C‘l{M). 

Note that RicM^ = has uniformly bounded norm, when measured 

via F^gk, hence when measured via g. Identities of the form (2.0.7)-(2.0.8) imply 
F^ RicMfc —^ RicM in for each £ > 0, p < oo. But the observation above 

implies some subsequence converges weak* in L°°. It follows that 

RicM e L“(M). 

A similar argument gives Ric^M ^ L°°{dM), and also similarly we obtain for the 
mean curvature H of dM ^ M that H G Lip{dM). Thus Theorem 2.1 applies 
to give g G C^{M), in boundary harmonic coordinates. This finishes the proof of 
Theorem 3.1. 

Remark. Invoking the definition of the Gromov-Hausdorff topology (cf. [Gr]) we 
can show that M(Ro, ig, Sg, dg) is compact in the Gromov-Hausdorff topology and 
C'’'-convergence is equivalent to Gromov-Hausdorff convergence on this compact set, 
for any r G [1, 2). 
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§3.4: Convergence of the geodesic flow and implications 

Here we establish results on the limiting behavior of geodesic hows under 
convergence, and implications for the injectivity radius, that sharpen and generalize 
some of the results of [Sak], We work in the following setting. Let M be a hxed 
compact manifold, with a C* coordinate system. Let g, gk be metric tensors on M. 
Assume 

(3.4.1) g,gkeCl(M), gk ^ g in C\M). 

Say these metric tensors dehne exponential maps 

(3.4.2) Expp : TpM D U ^ M, Exp^i, ^ : TpM D U ^ M, 

where U is a neighborhood of 0 G TpM if p G and U consists of n G TpM in 

a neighborhood of 0 that point into transversahy to dM, if p G dM. Here is 

our hrst result. 

Proposition 3.4.1. Under the hypotheses listed above, given v E U C TpM , 

(3.4.3) lim Exp;t,p('^) = Exp (n). 

k^oo 

Proof. Let Xk,X denote the vector helds on TM generating these hows, essentially 
the Hamiltonian vector helds associated with gk and g. Thus the coefficients of 
these vector helds belong to Cl and there is C^-convergence Xk —^ X. Orbits of 
Xk are 

(3.4.4) Pfc(t) = ( 7 fc(t), 7 [(t)), 

and we have p(.(t) = Xkigjkif)) uniformly bounded, hence both 7 ^ and 7 j(. uniformly 
Lipschitz. It follows from Arzela’s theorem that a subsequence converges locally 
uniformly: 

(3.4.5) 7 fc^ a, 7 fc_^ a'. 

Clearly (a, a') is an orbit of X, and a(0) = p, cr'(O) = v. Since (by Osgood’s 
theorem) the how generated by X is unique, it follows that a = 7 . Since this is 
true for any convergent subsequence, the result (3.4.3) follows. 

We next discuss the injectivity radius. We need to modify the dehnition used for 
(3.0.2), since in the present case Exp^ : U ^ M need not be or even Lipschitz. 
If p G and p < dist(p, dM), we will say 

(3.4.6) i{P:g)>P 
provided that, for each unit vector v G TpM"^"^^, 

(3.4.7) |t| < p ^ dist(p, Expp(tu)) = |t|, 

i.e., provided 7 t;(s) = Expp(su) is length-minimizing from p to 'jvit) for |t| < p. 
This does imply injectivity: 
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Proposition 3.4.2. If i{p, g) > p, then 

(3.4.8) Expp : i?p(0) —M*”* is one-to-one^ 

where i?p(0) = G : g{v,v) < p^}. 

Proof. If f, w G TpM^'^^ are unit vectors and 7 i;(s) = Jwit) with |s|, \t\ < p, then the 
condition (3.4.7) forces s = t (maybe after changing the sign of rc). If n ^ w, these 
geodesics must intersect non-tangentially aX q = 'jvit) = he., at a positive 

angle. Then a standard construction produces, for small £ > 0, a curve from p to 
7 t,(t+ £) shorter than \t\ + £, contradicting (3.4.7). 

We next have the following semi continuity result. 

Proposition 3.4.3. In the setting of Proposition S.f.l, 

(3.4.9) f(p, pfc) > po V /c f(p, g) > po. 

Proof. Given v G TpM'^'^^ such that g{v,v) = 1, we set Vk — v/^/gk{v^, so 
Qkivk, Vk) = 1 and Vk —> v. We are given that for all /c, 

(3.4.10) \t\ < Po ^ distfc(p, Expfc p(tnfc)) = \t\, 

where distfc denotes distance as determined by the metric tensor gk. Now it is 
elementary that 

(3.4.11) distfc(p, q) —^ dist(p, q), as /c —> cx), 
uniformly for q E M. Hence 

(3.4.12) distA;(p, Expp(tn)) ^ dist(p, Expp(fn)), 
while Proposition 3.4.1 together with (3.4.11) implies 

(3.4.13) distfc(p, Expfc p(tnA;)) - distfc(p, Expp(tn)) ^ 0. 
so we have 

(3.4.14) \t\ < Po dist(p, Expp(fn)) = \t\, 
as desired. 

We can also dehne the following sort of “boundary injectivity radius.” We say 

(3.4.15) lt(M,g)>p 

provided that for each p G 9M, inward normal Up G TpM (orthogonal to Tp(dM)), 

(3.4.16) 0 <t < p dist(Expp(fr'p), dM) = t. 

Parallel to Proposition 3.4.2, we have: 
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Proposition 3.4.4. > p, then 

(3.4.17) $:aMx[0,p)—^ M, 
given by 

(3.4.18) $(p, t) = Expp(tz/p), 
is one-to-one. 

Then, parallel to Proposition 3.4.3, we have 
Proposition 3.4.5. In the setting of Proposition S.f.l, 

(3.4.19) ih{M, gk) > po V /c ib{M, g) > po- 

The proofs of these resnlts are very similar to those of their connterparts above. 
Finally, we have the following signihcant implication for the geometric convergence 
obtained in Theorem 3.1. 

Corollary 3.4.6. If{M,g) is a limit of {Mk, gk) ^ M(Ro, io,‘So, do) as in Theorem 
3.1, then 

i{p, g) > niin(io, dist(p, dM)), V p G 

(3.4.20) i{p, g\dM) > *o, V p G dM, 

k{M,g) > io- 


4. Gel’fand inverse boundary problem 

In this section we prove nniqneness and stability and provide a reconstrnction 
procednre for the inverse bonndary spectral problem. To hx notations, assnme that 
(M, g, dM) is a compact, connected manifold, with nonempty bonndary, provided 
with a metric tensor g with some limited smoothness (specihed more precisely be¬ 
low). Let be the Nenmann Laplacian. Denote by 

0 = Ai < A 2 A ■ ■ ■ 

the eigenvalnes (connting mnltiplicity) of —and by 

(/>i=Vo1(M)-i/2, 

the corresponding orthonormalized eigenfnnctions. 
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The Gel’fand inverse boundary problem is the problem of the reconstruction of 
{M,g) from its boundary spectral data, i.e., the collection (9M, {Xk, <^k\dM}’^i)- 

Remark. The data used in the original formulation of the Gel’fand inverse bound¬ 
ary problem [Ge] consists of the trace on dM of the resolvent kernel, R\{x, y) of 
given for all A G C\spec A^, x, y G dM. The equivalence of these data and bound¬ 
ary spectral data, and also dynamic inverse boundary data for the corresponding 
wave, heat and non-stationary Schrodinger equations is proven in [KKL], [KLM], 
as are analogous results for the Dirichlet problem. 

The main results of this section are the uniqueness result below and stability 
result in §4.3. 

Theorem 4.1. Let M be a compact, connected manifold with nonempty bound¬ 
ary and metric tensor. Then the boundary spectral data {dM, {Xk, (t>k\dM}'^i) 
determine the manifold M and its metric g uniquely. 

Such a result was established in the case in [BKl], taking into account 
[Ta]; see also [KKL]. Our proof here will incorporate techniques from these papers, 
plus some additional arguments necessary to handle the reduced smoothness. One 
essential role played by the G^-hypothesis is that this implies non-branching of 
geodesics on M, including geodesics passing transversally from dM. In §4.1 we 
show that M is uniquely determined as a topological space. We proceed in §4.2 to 
show that the differential structure and metric tensor on M are uniquely determined. 
In §4.3 we apply Theorem 4.1 together with Theorem 3.1 to establish a result on 
the conditional stability of this inverse boundary problem, building on results of 
[K2L]. 


§4.1: Determining the domain 

We start with the introduction of some useful geometric objects. Let T C dM 
be open and take t > 0. Then we set 

(4.1.1) M(r, t)^{xeM : d{x, T) < t}, 

the domain of influence of T at “time” t. This terminology refers to the correspond¬ 
ing wave equation where M{r,t) is the subdomain of M hlled by the time t with 
waves sent from T. 

Now let r consist of a hnite number Ti,..., T^, of subsets T and t"*", t~ be two 
m-dimensional vectors with positive entries, tf~ = {t ^,..., tff), t~ = {tf ,..., ). 

Then set 


M(r, t+, r) = []{M{r„t+)\M{r„t-)) c m, 

i=l 


(4,1.2) 
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and define 

(4.1.3) L(r, t+, D = j^L2(M(r, t+, D) c 

Here, T stands for the Fonrier transform of fnnctions from LP‘{M), 


(4.1.4) J^('U) = e U{x) = ^Uk(l)k{x), Wfc = (W, (/)fc)L2(M), 

k=l 

and the snbspace L^(M(r, t~)) consists of all fnnctions in LP‘{M) with snpport 

in the set M(r, t~). 

Two basic ingredients for the reconstrnction of the manifold M are the approxi¬ 
mate controllability and Blagovestchenskii’s formnla. 

The controllability result is an implication of Tataru’s unique continuation result 
for the wave equation ([Ta], see also [Ho], [Ta2]). To describe it, consider the wave 
equation 


{d^ — A)u-^{x, t) = 0 in M X M_|_ 

(4.1.5) f f f ^ 

w^|t=o = o, 'u/|t=o = o, iVw^|9MxR+= / e (^(((r X (o,T)), 

where N is the exterior unit normal held to dM. Using Tataru’s theorem, it was 
shown in Theorem 3.10 of [KKL] that the following holds. 

Proposition 4.1.1. For each T > 0, the set {u-^{T) : f G L^(r x (0,T))} is a 
dense subspace o/L^(M(r, T)). 

(Actually Theorem 3.10 of [KKL] is written to address the Dirichlet boundary 
condition, but the same argument works for the Neumann boundary condition.) 

Blagovestchenskii’s formula gives the Fourier coefRcients u(,{t) of a wave u-^{■ ,t) 
in terms of the boundary spectral data, 

(4.1.6) = [ [ /(a^UO —— —4>k{x)dSgdt'. 

Jo J V J^k 

dM 

To prove this one starts with df {u(t) , (f)k)= {Au{t), (f)k)L^ and applies Green’s 
formula to get an inhomogeneous ODE for {u{t)^4>k)L^^ yielding (4.1.6). 

Note that in the formula (4.1.6) there appears the Riemannian volume dSg of 
dM, which we are not given. However, we are given dM as a manifold, so an 
arbitrarily chosen volume element element has the form dS = ndSg, where k is 
G^-smooth and strictly positive. We can construct the Fourier coefficients of the 
wave u^^{x,t) for any boundary source /, and despite our lack of knowledge of k, 
we do have the following. 
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Corollary 4.1.2. Given F C dM and t > 0, the boundary spectral data determine 
the subspace 

(4.1.7) L(r, t) = t)) C f. 

In fact, let {f^, : v G } have dense linear span in L^(r x (0,T)). Then L(r,t) is 
the closed linear span in P' of '■ ^ ^ Z}? where G is given by pi,^k = '^k' if)- 
Thus we can find the orthoprojection P : P ^ L(r,t) to this subspace. 

From here, using the elementary identities 

(4.1.8) =P|L2(5,), L\A,\B,)=L\A,)nL‘^{B,)^, 

i i 

we deduce that the boundary spectral data uniquely determine the subspaces L(r, 
t~) of for any F, t~ with arbitrary m (compare [KKL] and [Bel]). In par¬ 
ticular, for any F, t+, t~ we can see if L(r, t+, t~) = {0} or not. Equivalently, we 
can see if M(r, t"*", t~) contains an open ball or not. 

Next, let h G C{dM). We can ask if h is the boundary distance function for 
some X G M. To this end, we choose points Zj G dM, j = 1,..., m, their small 
neighborhoods Fi,..., F^ and numbers = h{zj) ± 1/m. When m ^ oo, the fact 
that L(r, t~) 7 ^ {0} for any m determines whether there is a point a: G M such 
that h{z) = dist( 2 ;, x), z G dM. Thus we have shown that the boundary spectral 
data determine the image in L°°{dM) of the boundary distance representation R. 
Here, R : M ^ C{dM) is dehned by 

(4.1.9) R{x)=rx{-), ra;(z) = dist(a:, z), zEdM. 

(Compare [KKL] and [Kul]). Clearly, the map R is Lipschitz continuous. Moreover, 
under the assumptions of Theorem 4.1 it is injective. To see this, let and 

let z G dM be a point of minimum of these functions. Then both x and y lie on the 
normal geodesic to dM starting in ^ at the same arclength rx{z) = ry{z). As the 
metric g G C‘l{M), it follows from Corollary 2.5.2 that this normal geodesic does 
not branch. Therefore, x — y. 

Since M is compact, injectivity and continuity imply that i? is a homeomorphism, 
i.e., R{M) with the distance inherited from L°°{dM) and (M, ^f) are homeomorphic, 
and thus R{M) can be identihed with M as a topological manifold. We have 
established the following. 

Proposition 4.1.3. Assume {Mi,gi) and (M 2 , ( 72 ) satisfy the hypotheses of Theo¬ 
rem 4 . 1 . If they have identical boundary spectral data, including dMi = dM 2 = X, 
as manifolds, then there is a natural correspondence of R{Mi) and R{M 2 ) C 
C{X), producing a uniquely defined homeomorphism 

M 2 . 


(4.1.10) 


X : Ml 
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4.2: Determining the metric 


Our next goal is to reconstruct the differential and Riemannian structures on 
M. To this end, let us return to L(r, t~) and consider the orthoprojection 
P(r, t"*", t~) of £2 onto this subspace. Then, 


(4.2,1) 


(r(r, f+, { = 




(f)i{x)(f)j{x) dVx, 


where Cj = (0,..., 0,1, 0,...) is the sequence having a 1 at the jth place. Also, 


(4.2.2) 




(j)i{xfdV^ 


voi(M(r, t+, r)) 

Vol(M) 


Next choose a sequence (T^i., /c = 1, 2,..., with nik cxd, where is 

the dimension of so that M(r;., shrinks to {x} when k ^ 00 . Then by 

formulae (4.2.1)-(4.2.2) we see that 


(4.2.3) 


lim {P{T,^, t+, 4)ei, ej) ■ {P{T^, t+, tjei, d) 


- 1/2 


(f)j{x) 


Thus we can find values of the eigenfunctions (f)k{x) for all /c = 1, 2,... and x G M. 

To proceed further we need an auxiliary statement about the properties of the 
eigenfunctions. Let $ be the space of all finite linear combinations of (pk, k = 

1 , 2 ,.... 

Lemma 4.2.1. Under the assumptions of Theorem f.l, $ is dense in the space 
{u e H^’P{M) : Nu\dM = 0} for any s G [0,3), p G (1 ,cxd). Moreover, if x ^ 
there are n indices k{l ),..., k{n) (depending on x) and a neighborhood U of x such 
that (pk{i){x), ..., (pk{n){x) form a C^-smooth coordinate system in U. 

Proof. Assuming that g G , r G (1, 2), consider, for any p G (1, cxd), the Neumann 
Laplacian, A^, with domain 

(4.2.4) V{A^) = {ue : Nu\dM = 0}. 

Denote by t > 0 the corresponding contraction semigroup and by (—A^)^, s > 

0, the real powers of — A^, defined for s G (0,1) via subordination. By Stein’s 
Littlewood-Paley theory for symmetric diffusion semigroups (cf. [St]), 

(4.2.5) = P((-A^)^/2), s g [0, 00 ), 

is a complex interpolation scale, in s, for each p G (1, cxd). In particular, 

(4.2.6) T>((-A^)"/2) = H^’P{M), for 0 < s < 1. 
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Hence, for 0 < s < 1, 

^ P((-A«)‘+'/2) = {u e V(-A^) : An 6 ®((A~)‘/^)} 

= {ue : Nu\dM = 0, and An G H^’P{M)}. 

If we reqnire g G C^{M), we can nse regnlarity resnlts to obtain that 

(4.2.8) V{{-A^y+^/^) ^ {ue : Nu\dM =0}, 0 < s < 1. 

Since (1 — Lp{M) is an isomorphism, acting bijectively on $, 

the desired density of $ will follow from the density of $ in Lp{M). We now 
demonstrate this density. 

Snppose / G L‘^{M) (with q = p') and (/,«.) = 0 for all n G $. If Q > 2 then 
/ G LF‘{M) and clearly / = 0. So we need only worry abont the case q < 2. Note 
that (/, n) = 0 implies 

(4.2.9) u) = 0, V t > 0, n G $. 

Now is a holomorphic semigronp on so for all t > 0, 

(4.2.10) / G V{A^) C C 

with q 2 > q, by the Sobolev embedding theorem. Iterating this and nsing the 
semigronp property gives 

(4.2.11) / G L'?(M) ^ G L2(M) V t > 0. 

Hence (4.2.9) implies = 0, for all t > 0. Bnt ^ / in L^{M) as t \ 0, 

so / = 0. This completes the proof of the hrst statement of the lemma. 

To demonstrate the second statement, we hrst note that since s < 3 and p < oo 
are arbitrary, it follows from the hrst part of the lemma that 

(4.2.12) C closnre of $ in C^(M). 

Let now x G and (x^, • • • ,a:"^) be some local coordinates near x. Denote by 

Ta; : $ —the map, 

(4.2.13) T^iu) = idiu{x), ■ ■ ■ , dnu{x)). 

It follows from (4.2.12) that Ta;($) = M"-, i.e., there are indices /c(l),..., k{n) (de¬ 
pending on x) snch that V4>k{i) (3^)? * = 1, • • •, n, are linearly independent. Moreover, 
we know that the eigenfnnctions (f)k G This proves the second statement 

of the lemma. 

Having this, we are in a position to rehne onr statement abont the homeomor- 
phism X- Ml ^ M 2 established in Proposition 4.1.3. 
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Proposition 4.2.2. Let {Mi,gi) and (M 2 , ( 72 ) be as in Proposition 4-P3, with 
identical boundary spectral data. Then the map x (4-1-10) has the property that 

(4.2.14) X ■ M^""* —M^""* is a C"^-diffeomorphism, and x *92 = 9i- 

Proof. We use the fact that, if {4>j} are the normalized eigenfunctions of on 

Ml and {<f>j} those on M 2 , then, as a consequence of (4.2.3), 

(4.2.15) ^^j{x{x)). 

Given p G ^ there exist indices /c(l),..., k{n) such that <f>k{i)^ • • ■, 4>k{n) form 
a local coordinate system on a neighborhood oi p = x{p)- Then, for x near p, we 
have 

(4.2.16) X ^ ((/)fc(i)(a:),...,(/)fc(n)(x)) = {^k{i){x{x)),.. . Ak{n){x{x))) ^ x(^) 

a composition of smooth maps. Thus x in (4.2.14) is smooth on a neighbor¬ 
hood of each p G hence on Ml^^. Interchanging the roles of Mi and M 2 , we 

have the same result for x~^- 

Finally, we show that the metric tensor is uniquely determined. For notational 
simplicity, just consider M = Mi. Let {x ^,..., x"^) be a coordinate system in a 
domain U C e.g., the one obtained earlier from the eigenfunctions. Then, for 

all /c = 1 , 2 ,..., 

(4.2.17) {x)didj(t)k{x) - P{x)di(t)k{x) = XkCpkix), P = 

where all eigenfunctions (fk and, henceforth, their derivatives as well as Ajt are 
already found. Let us consider equations (4.2.17) as linear equations for gP{x) — 
g^^{x), P{x). Using again (4.2.12), we see that the map ]R"'+^(^+i)/ 2 ^ 

(4.2.18) Tx{u) = {diu{x), didju{x) : i < j = 1, ■ ■ ■ ,m), 

is surjective (compare with (4.2.13)). Thus equations (4.2.17) are uniquely solvable 
since the 2 -jets of the eigenfunctions (pk at x span the whole space K"-+"-("-+i)/ 2 _ 

It follows that the diffeomorphism x ia (4.2.14) pulls the metric tensor 772 back 
to gi. The proof of Proposition 4.2.2, and hence of Theorem 4.1, is complete. 

We can get some more insight into how the geometry of (M, g) is determined by 
the boundary spectral data, particularly through (4.2.3), by examining further the 
maps 

(4.2.19) Tfc : M ^k{x) = ((/>i(a:),..., 

For simplicity we assume the eigenfunctions are arranged to be real valued. The 
argument proving Proposition 4.2.3 shows that for each compact K C M*”*, there 
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exists k such that \l/fc restricted to K is an embedding. In fact, we can do better 
than that, though for no k will \l/fc be an embedding of M, since D'^k{x) annihilates 
the normal to dM for each x G 9M, each k. Note that Lemma 4.2.1 implies the 
space of restrictions of elements of $ to dM is dense in C‘^{dM), so we can hnd 
ko = kQ{M, g) such that 

(4.2.20) \l/fc : dM is an embedding, for k > ko- 

We now augment 4/^ to 

(4.2.21) vl/f : M ^ R^+i, 4/#(x) = • • ■, 

where ^/^o G C^{M) is the eigenfunction for the Dirichlet problem, with smallest 
eigenvalue: 

(4.2.22) Ai^o = MdM = 

normalized by 

(4.2.23) J \Mx)\‘^dV = l, 

M 

and let us insist > 0 on Hopf’s principle implies N^/;Q{x) 0, V a: G 

dM. Hence the map (4.2.21) has the property 

(4.2.24) D'^f{x) is injective, y x e dM, k > ko. 

Thus injectivity holds on a collar neighborhood of dM in M. Combined with our 
previous observation about 4/^ embedding compact K C M*”*, this implies that 
there exists ki = ki{M,g) such that 

(4.2.25) D'^f{x) is injective, V a: G M, k > ki. 

Also we have 4/^ embedding dM. Perhaps increasing ki to be sure points are 
completely separated, we have the following result. 

Proposition 4.2.3. Let M he as in Theorem 4-2. Then there exists /c 2 = k 2 {M,g) 
such that 

(4.2.26) 4/^ : M —is a C'^-embedding, \/ k > k 2 . 

Having this, we can improve Proposition 4.2.2, as follows. 
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Corollary 4.2.4. In the setting of Proposition 4-2.2, we have 

(4.2.27) X : Mi — M 2 is a -diffeomorphism. 

Proof. Apply the construction (4.2.19)-(4.2.26) to obtain embeddings 

(4.2.28) :M2 —^ 

for k sufficiently large. Here we take \1/^ as in (4.2.21), with (M, g) = (Mi, ^fi), and 
dehne ^ 

= {'fo{x),^i{x),... ,^k{x)) 

on (M 2 , ( 72 ) in the analogous fashion. Using Proposition 4.2.2, we see that x pulls 
back the volume element of (M 2 , ( 72 ) to that of (Mi, 771 ), and deduce that 

(4.2.29) %fo{x) ^'ipo{x{x)). 

In concert with (4.2.15), this gives 

(4.2.30) m* = ^*ox. 

Since both maps in (4.2.28) are C^-embeddings (onto the same range), the result 
(4.2.27) follows from the implicit function theorem. 


Remark 1. One can readily strengthen the regularity results on x, given above, to 
regularity of class C^. Details can be left to the reader. 

Remark 2. The reason we do not use in this section the construction of [KKL] 
to recover the differential and Riemannian structure of (M, 77 ) is the following. In 
[KKL] we use as coordinates some distance functions on M. However, it is well 
known (e.g., [DTK]) that the resulting coordinates, in principle, lose two orders of 
regularity, i.e., in our case are just Cl regular. Clearly, it is rather inconvenient to 
work with such coordinates. 

Remark 3. There is an analogue of Theorem 4.1 in the case of the Dirichlet bound¬ 
ary spectral data (9M, {A/j, 77'0i|om}^i)- The proof is actually simpler. However, 
it requires some modihcations because we can no longer hud Vol(M(r, t"*", t~)). 
Therefore, instead of 'fkix) we hud 

= xeM^^\ k>l. 

However, an analog of Lemma 4.2.1 is valid and we can construct (PTsmooth co¬ 
ordinates in the vicinity of any x G M*"^* using functions fk, k > 1. Moreover, the 
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map Tj;, given by (4.2.18) is surjective on the linear span of the functions ^k- Since 
these functions satisfy equations 

-g"^{x)didj^k{x) - ¥{x)di^k{x) = {Xk - Xi)^k{x), 

with the same metric tensor but different 6% we can use these equations to hnd 
g'^^{x), X G Having found the metric inside M we can return to the analogs 

of equations (4.2.1), (4.2.2) and reconstruct i^kix), x G 


§4.3: Stabilization of the inverse problem 

In this section we consider stabilization of inverse problems using geometric con¬ 
vergence results and apply them to the Gel’fand problem. The basic thrust of our 
argument provides an illustration of a general “stabilization principle for inverse 
problems,” which we can describe abstractly as follows. 

Suppose A1 is a collection of objects, one element M of which you want to identify 
via the observation of data 'E’(M), in some set B of observable data. Suppose M. 
and B have natural topologies, and the map 

V:M —yB 

has been shown to be continuous. Suppose the uniqueness problem has been solved, 
so you know this map is one-to-one. However, typically such a map does not have 
a continuous inverse. This is a standard situation in the study of inverse problems, 
giving rise to the phenomenon of ill-posedness. This problem is made more acute 
by the fact that what one measures is not exactly equal to T>(M), but only an 
approximation to it. 

The key to the stabilization, which is useful for a wide variety of inverse problems, 
requires an a priori knowledge that the object M one wants to identify actually 
belongs to a subset M.q of Al, and that furthermore its closure in Af, Afo? is 
compact. In that case the restriction 

V:Mq — ^B 

is automatically a homeomorphism of AIq onto its range in B. Thus, when trying 
to identify the desired object M, one minimizes some measure of the difference 
between the calculated data V{Mj) and the observed data, while constraining Mj 
to belong to Afo- 

Having set up the abstract stabilization principle, we show how it applies to the 
Gel’fand problem. 

To prepare for this discussion, let us set up some notation. Denote by Mx(C^) 
the set of compact, connected manifolds M with nonempty boundary A, endowed 
with a metric tensor in C^{M). Given {M,g) G Afx(C'*), set 

(4.3.1) 1D(M,(7) = {A,,(/>,|x}°^i, 
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the right side denoting the bonndary spectral data of {M,g). We have 
(4.3.2) V : Mx(C^) Sx, 

where jBx denotes the set of seqnences : j > 1 }, with juj G /uj + cxd , 

and ipj G L‘^{X), modnlo an eqnivalence relation, which can be described as follows. 
We say if = aj'il^j^x) for some aj G C,\aj\ = 1. More 

generally, if gko = ■ ■ ■ = gki , we allow 


fci 

(4.3.3) ajk'^kix), j = kor--,h, 

k=kQ 


for a nnitary I x I matrix (ctj/c), I = ki — ko + 1. The content of Theorem 4.1 is that 
the map (4.3.2) is one-to-one. 

There are natnral topologies one can pnt on the sets in (4.3.2). On A4x{C^) 
one has the topology of convergence, for any r G (1,2), dehned in §3. On Bx 
one has a topology described as follows. We describe when {g’j, '0^} — {gj, i/jj}, as 
z/ —> cx). First we reqnire gj for each j. Next, if gk is simple, i.e., different 

from fXk-i and we reqnire a^i^k '^j some a'^ G C, \a^\ = 1. 

More generally, if Hk has mnltiplicity i, say Hk = ■ ■ ■ = fJ^k+e-i, we reqnire that 
there exist nnitary £ x £ matrices {a’:^j)k<ij<k+£-i snch that 

k+e-i 

(4.3.4) in L\X). 

j=k 


Compare [K2L]. 

Given these topologies, it follows from standard techniqnes of pertnrbation theory 
(cf. [K]) that V is continnons in (4.3.2). 

Now the map (4.3.2) is by no means invertible. This is a standard sitnation 
enconntered in the stndy of inverse problems, giving rise to the phenomenon of 
ill-posedness. One wants to “stabilize” the inverse problem, showing that certain 
a priori hypotheses on the domain {M,g) pnt it in a snbset K C Mx(C'^) having 
the property that can be shown to act continnonsly on the image of K. The 
resnlts of §3 provide a tool to accomplish this. 

Recall the class M(Ro, ioj 5'o, do) dehned in §3. Given a bonndary X, let Mx(Roj 
io, So, do) denote the set of snch manifolds with bonndary X. It follows from The¬ 
orem 3.1 that A4x(Ro, k, So, do) is compact in the topology, for any r G (1, 2), 
and is contained in M.x{C‘l). We hence give M.x{Ro,io, So,do) the (7” topology, 
and we see this is independent of r, for r G (1, 2). 

Combined with Theorem 4.1, these observations yield the following conditional 
stability of the Gel’fand inverse problem. 
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Theorem 4.3.1. Given Rq, io, Sq, do G (0, cx)), 

V : A4x(Ro, io, ‘So, do) — ^ Rx 

is a homeomorphism of Aix{Ro, io, So, do) onto its range, Bx{Ro, io, So, do); hence 
V ^ : BxiRoi ioi So, do) —^ Mx(Ro, *o, So, do) 


is continuous. 

Thus, if {Mk,gk), iM,g) G fAxiRoGo, So,do) and the boundary spectral data 
of {Mk,gk) tend to the boundary spectral data of {M,g) in B, then, for large k, 
Mk are diffeomorphic to M and gk ^ g , for all r < 2. 


5. Auxiliary regularity results 

Here we establish a number of elliptic regularity results, needed in the analysis 
in §2, which we did not hnd in the literature. These results tend to be variants 
of known results, but they differ in various key respects. Sometimes it is in the 
category of function space involved, e.g., coefRcients of the PDE in a non-standard 
space, which nevertheless arose naturally in the Ricci equation analysis. In some 
cases we can get away with a short argument based on standard results, while in 
other cases we need to do more work. 

In §5.1 we establish a local regularity result for an elliptic PDE with coehcients 
simultaneously satisfying a Holder condition and a Besov condition. In §5.2 we 
establish estimates for the Dirichlet problem when the coefficients simultaneously 
satisfy a Holder condition and a Sobolev condition. In §§5.3-5.5 we obtain estimates 
on weak solutions to Neumann boundary problems, with rough coefficients. 


§5.1: Local Besov regularity 

The following result establishes (2.2.8). 

Proposition 5.1.1. Assume u G H^’‘^{0) solves the elliptic PDE 

(5.1.1) djtt^^dkU = 0 on O. 

Assume 

(5.1.2) G C" n p 

with r, s E (0,1), p G (1, cxo). (One should assume r < s.) Then, locally, 

(5.1.3) u G n Bl+T 
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Proof. That u G is well known; we show u G ■ The proof is like that of 

Proposition 9.4 in Chapter III of [T2]. In particnlar, we nse paraprodncts, operators 
of the form TaU, a tool in nonlinear PDE introdnced by J.-M. Bony. A sketch of 
the behavior of paraprodncts can be fonnd in Chapter II of [T2]. 

Set = djT^jkdk and write (5.1.1) as 

(5.1.4) L*u = -dj[Td^ua^^ + , d^u)]. 

Here G OPBSf i is elliptic and, since G C'’’, by Proposition 6.1 in Chapter I 
of [T2] we have E G OPS^l snch that EL^ = / + F, F G OPS^l- Then we have 

(5.1.5) u = -Edj[Td^ua^^ + dku)] - Eu. 

Now 

u G (7^+" ^ Td,u: Rduu e OPSl^ 

(5.1.6) ^ Td,ua^^ + F(a^■^ dku) G 

^ Fc>,[T9,„a^" + F(a^", dku)] G F;+^ 

the second implication nsing the hypothesis that G B^ ^. Thns (5.1.5) gives 

(5.1.7) u = —Eu mod 

and since F G OPSff \ an iteration from u G readily yields u G B^'^^. 

§5.2: FP-SoBOLEV estimates for the Dirichlet problem 

Let O be a smooth, n-dimensional, manifold with bonndary, with metric tensor 

(5.2.1) g^k e C\n) n Fi’^'(O), s G (0,1), p G (1, oo). 

We aim to prove the following regnlarity resnlt. 

Proposition 5.2.1. Let O C dPt he open. Assume u G Lip(O) satisfies 

(5.2.2) Au = f e LP{n), u\^ = 0. 

Also assume p>n/{n — l). Given z E O, there exists a neighborhood U of z in Pi 
such that 


( 5 . 2 . 3 ) 


u G H^'P{U). 
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Proof. Take ip G such that p = 1 near ^ and </; = 0 on a neighborhood in O 

of dfl \ O. Suppose p is supported in a coordinate patch. With , we 

have 


(5.2.4) dja^’^dkipu) = pg^^'^f + u{dja^^dkp) + 2a^'^{dkp){dju), 

and if (5.2.1) holds and u G Lip(O), each term on the right side of (5.2.4) belongs to 
L^iVt). This reduces us to the case fl = R!):, with = 0 and u having compact 
support in fl, where we relabel pu as u. (This part of the argument still works even 
if we weaken the hypothesis u G Lip(O) to n G L°°(f2) fl 

Computing formally, we have for 1 < £ < n — 1 that U£ = diu satishes 

(5.2.5) djtt^’^dkUi = di{g^l‘^f) - dj {{d(,a^^)dku ). 

We have g^^'^f G L'^{Q) and {d(,a^^){dku) G Lp( 0), under our hypotheses, so the 
right side of (5.2.5) belongs to i7“^’P(0), and = 0. 

We claim this implies 

(5.2.6) G if^’P(O), l<£<n-l. 


This follows from Theorem 5.5.5’ of [Mo2], as long as p > n/(n — 1). Granted this, 
we have djdkU G Lp{Q) for all j, k except j = k = n, and the standard trick of using 
the PDE (5.2.2) to solve for dfu yields dfu G Lp{Q), completing the proof. 

Noting some alternative conditions that imply the right sides of (5.2.4) and (5.2.5) 
belong to Lp{Q) and respectively, we have the following extension of 

Proposition 5.2.1. 


Proposition 5.2.2. The conclusion (5.2.3) holds for a solution to (5.2.2) provided 

(5.2.7) G ^"(n) n77^’“(0), ue 

with s G (0,1), a, 6 G (1, cx)), and 


(5.2.8) 


111 
a 0 p 


Returning to the setting of Proposition 5.2.1, we note the following simple corol¬ 
lary, which is directly applicable to establish (2.2.9). 

Corollary 5.2.3. The conclusion of Proposition 5.2.1 remains valid if (5.2.2) is 
generalized to 

(5.2.9) An = /GLP(0), u\^ = g E nUpiO). 

Proof. After perhaps shrinking O, we can assume g = G\o with G G i7^’^(0) fl 
Lip(O). Then u — G solves 

(5.2.10) A{u-G) = f E LP{n), n-G'|^=0, 
and Proposition 5.2.1 applies to n — G. 
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5.3: Regularity for weak solutions to the Neumann problem 


Let M be a smooth, compact, connected manifold, of dimension n. Assnme M 
has a Riemannian metric tensor that is Holder continnons, of class (7^, for some 
r G (0,1). Let O C M be a connected open set, with bonndary of class 
Actnally we can assnme dVl is smooth, since a diffeomorphism can smooth 

ont 90 while prodncing a new metric tensor of class C^. Let A denote the Laplace 
operator on M. Assnme V G R > 0 on M, and V > 0 on a set of positive 

measnre. Consider L = A — V. 

A weak solntion to the Nenmann problem 

(5.3.1) Lu = f on O, Nu = g on 90 

is an element u G iL^(O) satisfying 


(5.3.2) 


{du, d^jj) dV = — {Vu + f)'ijjdV — / g'lpdS^ 


n 


a 


9f2 


for all G (O). Here the volnme element dV on O and the area element dS on 
90 are determined by the Riemannian metric tensor on M, as is the inner prodnct 
(^,? 7 ) of 1-forms. We aim to prove the following: 

Theorem 5.3.1. Given s G (0, r), u G iL^(O) satisfying (5.3.2), p > n/(l — s), 
(5.3.3) / G LP(0), g G (7"(90) (7^+"(0). 


Onr hrst rednction is to show that it snfhces to take / = 0. Indeed, extending / 
by 0 on M \ O we can solve Lv = / on M. By Proposition 2.3 of [MT], we have 
V G (7^"''^(0), provided Lp{M) C which holds as long as C 

i.e., as long as p(l — s) > n. Then T|n satishes 

(5.3.4) Lv = f on O, Nv — go E (7^(90), 

and it snfhces to show that w = u — v belongs to (7^"*'^(0). 

Onr next step is to look at 

(5.3.5) Lw = 0 on O, Nw — gi E C^{dLl), 

where gi = g — go, and prodnce a solntion w E (7^"*'^(n). (If R = 0 on O, assnme 
JdQdi dS = 0.) Prodncing snch a solntion to (5.3.5) will prove Theorm 5.3.1, since 
a solntion w E iL^(O) is nniqne (np to an additive constant if R = 0 on 11). To 
prodnce snch a solntion, we nse the method of layer potentials. 
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Thus let E{x,y) be the integral kernel of L ^ : H ^(M) —> and dehne 

the single layer potential 

(5.3.6) Sh{x) = j E{x,y)h{y)dS{y). 

an 

To proceed we need an analysis of E{x,y). It is elementary to show that h G 
L^{dfl) Sh G \ 90). In fact, given compact E C 90, 

(5.3.7) h G L^(90), supp/i cE^She Cl+J{M \ E). 

This permits us to use partition of unity arguments and localize the study of Sh 
to coordinate patches. We can choose local coordinates such that 90 is given by 
{x : x'^ = 0}. 

As in [MT], we can write 

(5.3.8) E{x,y)^/g{y) = eo{x -y,y) + ei{x,y), 
where (if n > 3) 

\-(^- 2)/2 

(5.3.9) eo{x-y,y) =Cny2^gjkiy){xj -yj){xk-yk)j 
As shown in Theorem 2.6 of [MT], we have, for each £ > 0, 

\ei{x, y)\ < Ce\x — 

(5.3.10) " ' ^ 

\V^eiix,y)\<Ce\x-y\-^^-^--+^l 

Also (2.67) of [MT] implies (for 0 < s < r) 

(5.3.11) \'^xei{xi,y) - Vxei{x2,y)\ < Cs\xi - X2\'' \xi - 
provided \xi — X 2 \ < (l/2)|a:i — y\- 

Given h supported in E, the intersection of 90 with a coordinate patch, we 
analyze Sh as a sum of two pieces, Sh = St^h + Sih^ where 

Soh{x)= / eo{x - y,y)g{y)~^^‘^h{y)dS{y), 

9n 

Sih{x) = j ei{x,y)g{y)~^/‘^h{y)dS{y). 
on 


(5.3.12) 



Lemma 5.3.2. Given s G (0,r), we have 


(5.3.13) : L“(E) —^ C^+\M). 

Proof. We need to show that, for Xj G M, h G L°°(90), supported in E, 

(5.3.14) |ViSi/i(a:i) — VSih{x 2 )\ < C\xi — X 2 f ||h||Loo. 

There are two cases to consider. 

Case I. \xi — X 2 \ < (l/2)dist(a:i, E). 

Then use (5.3.11) to get (5.3.14). 

Case II. \xi — X 2 \ > (l/2)dist(a:i,E). 

Set C = {y G E : \xi—y\ < 4\xi—X2\}. Use (5.3.10) for y G C to analyze separately 
Va; Jq eo{xj — y, y)g{y)~^^‘^h{y) dS{y), and use (5.3.11) for y G E \ C to complete 
the analysis of the left side of (5.3.14). 

Lemma 5.3.3. Given s G (0,r), we have 

(5.3.15) 5o : C^{dn) —^ c^+^(n). 


Proof. This is a standard layer potential estimate. One has 

(5.3.16) IIV5o/i||Loo(f2) < C\\h\\c^{dn): 
and 

(5.3.17) \d‘^Soh{x',Xn)\ < C|a:n|^“^ ||/i||c-(an), 
which gives 

(5.3.18) \VSoh{xi)-VSoh{x2)\ < C\xi - X2\'' \\h\\cs(da), XjEVt] 

compare Proposition 8.7 in Chapter 13 of [Tl]. 

The next step in solving (5.3.5) in the form w = Sh is to analyze NSh\QCL- We 
have the standard formula 

(5.3.19) NSh\Q^={-h + K*y, 
with 


(5.3.20) 


K*h{x) = PV / N^E{x,y)h{y)dS{y), x G dQ. 
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Compare (2.81) of [MT] for such a formula in the more general context of a Lipschitz 
boundary dfl. We can break K* into two pieces, K* = Kq + Kl, using (5.3.8): 


(5.3.21) 


K*h{x) 


PV 


N^eoix 


an 


y,y)9iy) ^^‘^Hy)dS{y), 


and 

(5.3.22) 


Klh{x) 


j N^ei{x,y)g{y) ^/‘^h{y) dS{y). 

an 


By Lemma 5.3.2 we have : L°°(90) —> (7^(911). As for ATq, it is a pseudodif¬ 
ferential operator with double symbol, of the sort studied in Chapter I, §9 of [T2], 
with symbol 

(5.3.23) a{x,y,i)eC^S%. 

Furthermore, the appearance of Nx in (5.3.21) produces the following important 
cancellation effect: 


(5.3.24) a(a:, x, ^) = 0. 

(This is part of what makes analysis on domains with (7^^^ boundary easier than 
analysis on Lipschitz domains.) Hence, by Proposition 9.15 of [T2], Chapter I, we 
have Kq : L°°(90) —(7'^ (90), for all s < r. In summary, 

(5.3.25) A:*:L“(90)—^ (7"(90), V s < r, 
and hence 

(5.3.26) K* : C^{dQ) —(7^(90) is compact, for all s E (0,r). 

Thus 

(5.3.27) — I + K* : (7^(90) —(7^(90) is Fredholm, of index 0. 

2 

Using this we can prove the following. 

Proposition 5.3.4. Let gi E (7^(90), 0<s<r. If V > 0 somewhere (i.e., on a 
set of positive measure) on O, then (5.3.5) has a unique solution w = Sh E (7^"'“®(0), 
where 


(5.3.28) 


+ = hEC^idn) 
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IfV = 0 on Q, then (5.3.5) has a solution w G if and only if dS = 0, 

and such w is unique up to an additive constant. 

Proof. This is a standard argument in layer potential theory. One shows that, 
if y > 0 somewhere on O, then —{1/2)1 + K* is injective on C^{dPl), hence, by 
(5.3.27), bijective. IfV^ = 0 on bl, then —(1/2)7+77* has a one-dimensional kernel in 
(7® (90), and the constant function 1 annihilates its range. Compare the treatment 
of Theorem 3.4 in [MT], carried out in the more general context of a Lipschitz 
boundary. Once one solves (5.3.28) for h G (7^(90), the fact that w = Sh solves 
(5.3.5) and belongs to (7^"*'^(0) follows from the previous analysis. 

With this result, the proof of Theorem 5.3.1 is complete. 


5.4: Local regularity for the Neumann problem 


It is useful to strengthen the global regularity result of §5.3 to a local regularity 
result for a weak solution to the Neumann problem. With M, O, and L as in §5.3, 
let O be an open subset of 90 and suppose w is a weak solution to 

(5.4.1) Lu = f on O, = 9- 

That is to say, we assume u G 77^(0) and that (5.3.2) holds for all t/> G 77^(0) such 
that vanishes on a neighborhood of 90 \ (7 in 90. We prove the following. 

Theorem 5.4.1. Let u be a weak solution to (h.f.l). Assume u G 77^’'^(0) with 
q >2 and either q > 1/r or q > p. (Recall gjk G (7^./ As in Theorem 5.3.1, assume 
0 < s < r < 1 and p> n/{l — s). Then 

(5.4.2) / G LP(0), g G (7"((7) (7^+"(0 U O). 


Proof. Take pq & O and pick ip G C^{M), equal to 1 on a neighborhood of po, such 
that + = 0 on a neighborhood of 90 \ (7 in M. Let v = +^1^, so n G 77^’'^(0). We 
seek to establish extra regularity of n. It is readily verihed that n is a (global) weak 
solution of 


(5.4.3) 
with 

(5.4.4) 


Lv — f on O, Nv = g on 90, 


f = pf + 2.{dp, du) + u/S.p, 
g = pg + u{Np)\g^. 


We have 
(5.4.5) 


/ G LP^'?(0), g G (7"(90) + (7"(90) ■ Sj-i/^(90). 
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Extend / by 0 on M \ O and set vi = L ^f. By (2.16) of [MT], we have 

(5.4.6) C ^ vi e HP+^’P{M), p = r - e. 

In fact, if q > p, we can nse Proposition 2.3 of [MT] as in the beginning of the 
proof of Theorem 5.3.1, and say vi G C'^"*'^(M). If q < p, then we can take p > q in 

(5.4.6) . 

Now V 2 — V — vi\q, is n weak solntion to 

(5.4.7) Lv 2 = 0 on O, Nv 2 = §2 = 9 — Nvi, 
and 

(5.4.8) Nvi e C^dQ) ■ 

if q < p, and in (7®(90) if q > p. Note that pp > q{r — e) can be assnmed to be 
> 1 by the hypothesis rg > 1 , if we take £ > 0 small enongh. At this point we have 
92 G with 92 > q- From this we can dednce 

(5.4.9) ^2 G 92 >q. 

The resnlt (5.4.9) is a special case of mnch stronger known resnlts; let ns sketch 
the proof. First, parallel to (5.3.26), we have 

(5.4.10) K* : L^{dfl) — L^{dfl) is compact, for all p G (1, cxd). 

Indeed, the compactness of Kq follows from Proposition 9.5 in Chapter I of [T2], to¬ 
gether with (5.3.23)-(5.3.24), and the compactness of follows from the estimates 
in (5.3.10). Having (5.4.10), we see that 

(5.4.11) —/ -f K* : L^(90) — L^{dQ) is Fredholm, of index 0, 

for each p G (1, oo). Then an argument as in the proof of Proposition 5.3.4 yields 

(5.4.12) V 2 =Sh 2 , h 2 = + K*y \2 e yHdn), 

np to an additive constant if F = 0 on O. Now we have the non-tangential maximal 
fnnction estimate 

(5.4.13) ll(V‘5^2)*||L92(an) < 1 ^ 92 II^ 2 Hlw( an), 1 < ^2 < cxd, 

which is stronger than (5.4.9). The estimate (5.4.13), in the setting of (5.4.6) for 
a C'^'-metric tensor, is proven in (2.77) of [MT], in the more general context of a 
Lipschitz domain fl. 

Now (5.4.6) and (5.4.9) together give 

(5.4.14) n = -I- ^2 G i7^’'^^(0), 92 > 9- 

Now we go back to u. After shrinking O to a smaller neighborhood of po, we can 
replace onr hypothesis u G i7^’'^(0) by n. G i7^’'^2(0). Iterating this argnment yields 
u G (O) with q < 92 < 93 <■■■ ■ After a hnite nnmber of iterations we reach 
a point where Theorem 5.3.1 is applicable to n = (pu, and Theorem 5.4.1 is proven. 
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§5.5: Neumann data in 

We produce more regularity results on weak solutions to 

(5.5.1) Lu = f on O, Nu = g on dfl, 

starting with the following. As before, assume gjk € for some r > 0. 
Proposition 5.5.1. Let u G he a weak solution to (5.5.1). Assume 

(5.5.2) / G L^(0), g G L^{dVL), p > n, 1 < s < oo. 

Then we have 

(5.5.3) (Vn)* G L"(aO), n G (O). 

Proof. Since p > n, we can solve Lv = f (with / extended by 0) on a neighborhood 
OofO, withn G C^+'^{0), a > 0. This reduces our consideration to the case / = 0 
in (5.5.2). Then u is given by 

(5.5.4) u = Shj 
with 

(5.5.5) (^-h + K*^h = g, heL%dn). 

The fact that 

(5.5.6) h G L^dQ) ^ (VSh)* G L^{dn) 

is established in the context of a Holder continuous metric tensor (and in the more 
general context of a Lipschitz domain) in [MT]. 

The last part of (5.5.3) follows from the mapping property 

(5.5.7) 5 : L^{dn) —^ 77bnV(^-i)(0). 

This is demonstrated, in the more general context of a Lipschitz domain, in [MT2]. 
We describe here the basic structure of the argument. We write S = iSq + iSi, as in 
(5.3.12). Harmonic analysis techniques are brought to bear to establish 

(5.5.8) So : L ^{ dn ) —^ 

(Such a result is easier for a smooth domain than for a Lipschitz domain.) As for 
iSi, we already have in (5.3.13) that 

(5.5.9) : L“(aO) —^ 

for some a > 0. Meanwhile the estimate (5.3.10) on Vxei{x,y) is more than ade¬ 
quate to give 

(5.5.10) : L\dn) —^ 
and then interpolation gives more than 

(5.5.11) : L%dn) —^ 77bn./(n-i)^Q)_ 

We now establish a useful local regularity result. 
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Proposition 5.5.2. Let O be an open subset of dLl and assume u is a weak solution 
to 

(5.5.12) Lu = f on Vt, = 9- 

Take p > n, s G (1, cx)), and assume 

(5.5.13) ueH^’P{n), /eLP(O), geL^iO). 

Then each p E O has a neighborhood U in Vt such that 

(5.5.14) u E 

Proof. As in the proof of Theorem 5.4.1, we consider v = (pu, which is a global 
weak solntion of (5.4.3), with f,g given by (5.4.4). Under onr cnrrent hypotheses 
we have 

(5.5.15) /eLP(O), gEL^idn), 

so Proposition 5.5.1 gives v E and Proposition 5.5.2 is proven. 


References 

[Al] G. Alessandrini, Stable determination of condnctivity by bonndary mea- 
snrements, Applic. Anal. 27 (1988), 153-172. 

[Anl] M. Anderson, Convergence and rigidity of manifolds nnder Ricci cnrvatnre 
bonnds. Invent. Math. 102 (1990), 429-445. 

[An2] M. Anderson, Bonndary regnlarity, nniqneness, and non-nniqneness for AH 
Einstein metrics on 4-manifolds, Preprint, 2001. 

[Bel] M. Belishev, An approach to mnltidimensional inverse problems (Rnssian), 
Dokl. Akad. Nank. SSSR, 297 (1987), 524-527. 

[BKl] M. Belishev and Y. Knrylev, To the reconstrnction of a Riemannian manifold 
via its spectral data (BC-method), Comm. PDE 17 (1992), 767-804. 

[Bz] Y. Berezanskii, The nniqneness theorem in the inverse problem of spectral 
analysis for the Schrodinger eqnation (Rnssian), Trndy Moskov. Matem. 
Obsch., 7 (1958), 1-62. 

[Ch] J. Cheeger, Finiteness theorems for Riemannian manifolds, Amer. J. Math. 
92 (1970), 61-75. 

[CG] J. Cheeger and D. Gromoll, The splitting theorem for manifolds with non¬ 
negative Ricci cnrvatnre, J. Diff. Geom. 6 (1971), 119-128. 

[CW] Y.-Z. Chen and L.-C. Wn, Second Order Elliptic Eqnations and Elliptic 
Systems, Transl. Math. Monogr., No. 174, AMS, Providence, R.I., 1998. 



60 


[DTK] D. De Turck and J. Kazdan, Some regularity theorems in Riemannian ge¬ 
ometry, Ann. Scient. Ecole Norm. Sup. Paris, 14 (1981), 249-260. 

[Ge] I. Gel’fand, Some aspects of functional analysis and algebra, Proc. IGM 1 
(1954), 253-277. 

[GT] D. Gilbarg and N. Trudinger, Elliptic Partial Differential Equations of Sec¬ 
ond Order, Springer-Verlag, New York, 1983. 

[GW] R. Greene and H. Wu, Lipschitz convergence of Riemannian manifolds. Pa¬ 
cific J. Math. 131 (1988), 119-141. 

[Gr] M. Gromov, Metric Structures for Riemannian and Non-Riemannian Spaces 
(with appendices by M. Katz, P. Pansu, and S. Semmes), Birkhauser, Boston, 
1999. 

[Ha] P. Hartman, On the local uniqueness of geodesics, Amer. J. Math. 72 
(1950), 723-730. 

[HH] E. Hebey and M. Herzlich, Harmonic coordinates, harmonic radius, and 
convergence of Riemannian manifolds. Rend, di Mathem. Ser. VH, 17 
(1997), 569-605. 

[Ho] L. Hormander, On the uniqueness of the Gauchy problem under partial an- 
alyticity assumptions. In: Geom. Optics and Related Topics (F. Golombini 
and N. Lerner, eds.), Birkhauser, Boston, 1997. 

[KKL] A. Katchalov, Y. Kurylev, and M. Lassas, Inverse Boundary Spectral Prob¬ 
lems, Ghapman-Hall/GRG Press, Boca Raton, 2001. 

[KLM] A. Katchalov, Y. Kurylev, M. Lassas, and N. Mandache, Equivalence of 
time-domain inverse problems and boundary spectral problems. Preprint, 
2002 . 

[K] T. Kato, Perturbation Theory for Linear Operators, Springer-Varlag, New 
York, 1966. 

[Ka] A. Katsuda, BG-method and stability of Gel’fand inverse spectral problem, 
Proc. Gonf. “Spectral and Scattering Theory,” RIMS, Kyoto, 2001, 24-35. 

[K2L] A. Katsuda, Y. Kurylev, and M. Lassas, Stability in Gelfand inverse bound¬ 
ary spectral problem. Preprint, 2001. 

[Kod] S. Kodani, Gonvergence theory for Riemannian manifolds with boundary, 
Gompos. Math. 75 (1990), 171-192. 

[Ku] Y. Kurylev, Multidimensional Gel’fand inverse boundary problem and bound¬ 
ary distance map. In: Inverse Problems Related to Geometry (H. Soga, ed.), 
1-15, Ibaraki Univ. Press, Mito, 1997. 

[KuL] Y. Kurylev and M. Lassas, The multidimensional Gel’fand inverse problem 
for non-self-adjoint operators. Inverse Problems, 13 (1997), 1495-1501. 

[Mi] G. Miranda, Partial Differential Equations of Elliptic Type, Springer-Verlag, 
New York, 1970. 

[MT] M. Mitrea and M. Taylor, Potential theory on Lipschitz domains in Rie¬ 
mannian manifolds: Holder continuous metric tensors, Gomm. PDE 27 
(2000), 1487-1536. 

[MT2] M. Mitrea and M. Taylor, Local and global regularity on Lipschitz domains. 



61 


Manuscript, in preparation. 

[Mol] C. Morrey, Second order elliptic systems of differential equations, pp. 101- 
159 in Ann. of Math. Studies ^33, Princeton Univ. Press, Princeton, N.J., 
1954. 

[Mo2] C. Morrey, Multiple Integrals in the Calculus of Variations, Springer-Verlag, 
New York, 1966. 

[NSU] A. Nachman, J. Sylvester and G. Uhlmann, An n—dimensional Borg-Levinson 
theorem. Comm. Math. Phys. 115 (1988), 595-605. 

[Nv] R. Novikov, A multidimensional inverse spectral problem for the equation 
—A'0 + (v(x) — Eu{x))ip — 0 (Russian), Funkt. Anal, i Priloz. 22 (1988), 
11 - 22 . 

[P] S. Peters, Convergence of Riemannian manifolds. Compos. Math. 62 (1987), 
3-16. 

[Pe] P. Petersen, Riemannian Geometry, Springer-Verlag, New York, 1998. 

[Sak] T. Sakai, On continuity of iniectivity radius function. Math. J. Okayama 
Univ. 25 (1983), 91-97. 

[StU] P. Stefanov and G. Uhlmann, Stability estimates for the hyperbolic Dirichlet- 
to-Neumann map in anisotropic media, J. Funct. Anal. 154 (1998), 330-357. 

[St] E. Stein, Topics in Harmonic Analysis Related to the Littlewood-Paley The¬ 
ory, Princeton Univ. Press, Princeton, N.J., 1970. 

[Ta] D. Tataru, Unique continuation for solutions to PDE’s: between Hormander’s 
theorem and Holmgren’s theorem. Comm. PDE 20 (1995), 855-884. 

[Ta2] D. Tataru, Unique continuation for solutions to PDEs with partially analytic 
coefficients, J. Math. Pures Appl. 78 (1999), 505-521. 

[Tl] M. Taylor, Partial Differential Equations, Vols. 1-3, Springer-Verlag, New 
York, 1996. 

[T2] M. Taylor, Tools for PDE, AMS, Providence, R.I., 2000. 

[Trl] H. Triebel, Interpolation Theory, Function Spaces, Differential Operators, 
DVW, Berlin, 1978. 

[Tr2] H. Triebel, Function Spaces, Birkhauser, Boston, 1983. 


Michael Anderson, Mathematics Department, State Univ. of New York, Stony 
Brook, N.Y., 11794, USA 

Atsushi Katsuda, Mathematics Department, Okayama Univ., Tsushima-naka, 
Okayama, 700-8530, Japan 

Yaroslav Klfrylev, Department of Math. Sciences, Loughborough Univ., Lough¬ 
borough, LEll 3TU, UK 

Matti Lassas, Rolf Nevanlinna Institute, Univ. of Helsinki, FIN-00014, Finland 



62 


Michael Taylor, Mathematics Department, Univ. of North Carolina, Chapel 
Hill, N.C., 27599, USA 



